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Abstract 

Free  electron  laser  (FEL)  theory  is  extended  to  explain  several  effects  associated  with  high 
gain  operation.  The  trapped-particle  instability  is  reviewed  for  short  pulse  FELs,  FEL 
oscillators,  and  FEL  amplifiers.  A  new  FEL  theory  exactly  includes  the  effect  of  an  arbitrary 
electron  distribution  function  and  an  explanation  is  given  for  optical  self -guiding. 
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Introduction 

The  proposed  research  was  partially  funded  and  the  following  reports  on  the  resulting 
scientific  progress.  We  have  continued  to  work  closely  with  the  experimental  free  electron  laser 
(PEL)  projects  in  Universities  and  National  Research  Labs  with  emphasis  on  the  high-power 
configurations  at  LANL  and  LLNL  as  well  as  the  short  wavelength  configurations  at  Stanford  and 
LANL.  One  aspect  of  the  work  has  been  to  characterize  the  FEL  trapped-particle  instability.  An 
especially  important  result  is  a  new  theoretical  procedure  solving  the  FEL  wave  and  electron 
equations  in  weak  optical  fields  for  an  arbitrary  electron  distribution  function. 

The  slowly-varying,  self-consistent  wave  and  electron  equations  of  motion  provide  the  basic 
theory.  The  complex  optical  field  strength  is  a  and  the  electron  phase  is  Q?,x).  The  optical 
wave  envelope  is  then  described  by  the  parabolic  wave  equation: 

»  * 

~^l+^  a(?,x)  =  -<j(i)exp(-i$(*,x))>  ,  (1) 

•  « 

where  the  dimensionless  current  density  is  j(?,t)  and  is  the  transverse  Laplacian.  FELs 
amplify  the  radiation  with  a  co-propagating  relativistic  electron  beam  traveling  along  the  axis  of  a 
long,  undulating  magnetic  field.  For  the  periodic  undulator,  we  can  describe  the  electron 
dynamics  at  a  site  in  the  beam  with  the  self-consistent  pendulum  equation: 

C(z+T)  =  v(2+T)  =  -^[a(z)exp(iC(z+t))  +  a*{z)exp(-i;(2+x))]  .  (2) 

The  slippage  of  the  optical  field  with  respect  to  the  electrons  is  indicated  by  ( z  +  x )  .  The 
longitudinal  variables  have  been  normalized  to  the  slippage  distance  and  all  variables  have  been 
reduced  to  a  meaningful  dimensionless  form. 

While  simple  in  form,  these  equations  contain  a  large  array  of  important  effects  describing 
FEL  operation.  The  goals  of  high  optical  power  or  short  optical  wavelengths  both  have  in  common 
the  need  for  increasing  the  current  density  j(f)  above.  This  coupling  between  the  electrons  and 
light  can  be  increased  by  increasing  the  undulator  length  and/or  the  physical  electron  current 
density. 
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Panel  Discussion  on  Long  Undulators 

Several  members  of  the  FEL  community  reviewed  the  possible  limitations  of  long  undulators 
containing  up  to  103  periods.  The  panel  members  were  K.  Halbach,  B.  Kincaid,  B.  Newman,  D.  A. 
Q.  Deacon,  and  myself  (panel  discussion  leader).  The  general  conclusion  was  that  the  expected 
improvements  in  the  design  of  currently  proposed  undulators  will  make  it  possible  to  reach  the 
long  lengths  desired.  The  tasks  will  be  difficult,  but  not  forbidden  by  serious  physical  limitations. 
The  reference  is 

"Long  Wigglers  for  Free  Electron  Lasers:  A  Maximum  Limit?”,  First  International  Laser 
Science  Conference,  Dallas  TX  (1985). 

FEL  Review  Article 

Some  work  on  a  major  review  article  fell  into  this  contracting  period.  FELs  were  reviewed 
starting  with  their  historical  development  from  microwave  tubes  and  atomic  lasers.  Data 
describing  all  operating  FEL  experiments  was  presented  in  a  table  format.  The  basic  FEL  theory 
is  reviewed  with  extensions  to  include  many  of  the  complicated  effects  in  equations  (1)  and  (2). 
The  high-gain  growth  rate  is  derived  for  the  first  time  from  the  single-particle  equations  of  motion. 
Specific  accelerators  suitable  for  driving  FELs  are  reviewed  in  a  table  format.  The  linac  driven 
FEL,  induction  linac  driven  FEL,  and  the  storage-ring  FEL  are  described  in  more  detail.  The 
extensions  of  the  FEL  undulator  designs  and  new  configurations  are  summarized.  The  reference 
is 

W.  B.  Colson  and  A.  M.  Sessler,  "Free  Electron  Lasers",  Annual  Reviews  of  Nuclear  and 
Particle  Science  35,  25  (1985). 

A  document  report  on  this  paper  follows. 

The  Trapped-Particle  Instability 

Several  simulations  and  simulation  techniques  have  been  developed  describing  the 
trapped-particle  instability  in  free-electron  laser  oscillators  and  amplifiers.  There  are  distinct 
differences  in  the  physics  in  these  two  cases.  The  electrons  in  a  high-power  free-electron  laser 
can  become  trapped  and  oscillate  in  deep  potential  wells  formed  by  the  combined  optical  and 
undulator  field  forces.  The  trapped  current  oscillates  at  the  synchrotron  frequency,  and  can  drive 
the  optical  wave  at  sideband  frequencies  around  the  fundamental.  This  "trapped-particle 
instability"  can  occur  in  both  the  oscillator  and  amplifier  configurations.  The  gain  in  the  trapped- 
particle  sideband  frequencies  has  been  calculated  at  saturation  and  agrees  well  with  the 
simulation  results.  General  features  of  the  trapped-particle  instability  in  free-electron  laser 


oscillators  and  amplifiers  have  been  summarized.  Dimensionless  parameters  are  used  to  clarify 
the  trends  for  a  wide  range  of  FEL  designs.  The  instability  was  first  observed  in  the  simulation  of 
short  pulse  oscillators,  and  then  extended  to  long  pulse  oscillators.  In  this  research,  the  trends  of 
the  trapped-particle  instability  in  short  pulse  oscillators,  long  pulse  oscillators,  and  now  long  pulse 
amplifiers  is  compared.  This  research  made  the  first  simulations  of  the  sideband  instability  in  the 
LLNL  long  pulse  amplifiers.  There  are  two  publications: 

"The  Trapped-Particle  Instability  in  Free-Electron  Laser  Oscillators  and  Amplifiers",  W.  B. 

Colson,  Proc.  of  the  1985  Free-Electron  Laser  Conference,  Lake  Tahoe  CA  (1986). 

"The  Effect  of  Electron  Trapping  in  Free-Electron  Laser  Oscillators  and  Amplifiers",  W.  B. 

Colson,  Proceedings  of  the  1985  International  Conf.  on  LASERS,  Las  Vegas  CA  (1986). 
These  two  documents  are  included  at  the  end  of  this  report. 

Electron  Beam  Energy  Spread  and  Emittance 

In  an  FEL,  maintaining  the  coherence  of  the  electron  bunches  over  a  significant  interaction 
length  imposes  important  restrictions  on  the  electron  beam  quality.  An  energy  or  angular  spread 
(due  to  emittance)  contributes  a  random  component  to  the  electron  motion  that  decreases  the 
coherent  bunching  in  time.  Some  of  the  earliest  FEL  experiments  used  electron  beams  that  were 
essentially  monoenergetic,  but  practically  all  subsequent  experiments  have  made  use  of  higher 
current  sources  with  significant  energy  spread  or  emittance.  Many  accelerators  present  a  design 
trade-off  between  high  current  and  high  beam  quality.  This  makes  it  essential  to  accurately 
evaluate  the  effects  of  beam  quality  in  present  and  future  experiments.  It  is  particularly  important 
for  FELs  designed  to  operate  at  XUV  or  X-ray  wavelengths. 

The  theory  presented  here  uses  a  convenient,  yet  powerful,  method  of  including  an  arbitrary 
electron  distribution  function  in  a  self-consistent  integral  equation  for  the  complex  optical  field. 
FEL  gain  and  the  effects  of  beam  quality  can  then  be  calculated  analytically  or  integrated  on  a 
small  computer.  Since  the  basic  equations  solved  here  are  the  same  as  in  computer  simulations 
or  the  plasma  dispersion  methods,  specific  physical  results  have  been  shown  to  agree  with  those 
methods  when  a  direct  comparison  is  possible.  The  computer  simulations  have  proved  to  be  a 
useful  method  of  understanding  many  aspects  of  the  FEL  interaction,  but  one  of  the  most  difficult 
effects  to  accurately  characterize  is  that  of  electron  beam  quality.  Even  a  prohibitively  large 
number  of  sample  particles  is  far  short  of  the  number  in  a  real  experiment,  and  yet  introduces  a 
large  amount  of  numerical  noise  when  distributed  over  a  large  volume  of  phase-space.  To 
reproduce  some  of  the  results  shown  later  in  this  paper,  we  found  the  simulation  method  to  be 
several  hundred  to  a  thousand  times  less  efficient.  While  many  other  FEL  topics  are  most 


efficiently  studied  through  simulations,  the  detrimental  effects  of  beam  quality  are  probably  better 
handled  through  a  combination  of  analytic  and  numerical  techniques. 

The  stability  analysis  used  to  obtain  plasma  dispersion  relations  usually  calculates  the 
reduced  FEL  growth  rates  due  to  poor  beam  quality.  This  method  can  lead  to  analytical 
expressions,  but  depends  upon  specific  models  for  the  electron  beam  distribution,  and  does  not 
easily  describe  more  complicated  transient  behavior  where  the  FEL  growth  rate  is  not  constant; 
the  FEL  is  often  designed  to  operate  in  this  regime.  In  addition,  the  exact  formulation  presented 
here  works  smoothly  between  different  regimes  of  operation  like  high  and  low  gain.  The  only 
requirement  is  weak  optical  fields. 

An  FEL  with  arbitrary  beam  quality  is  accurately  described  by  the  Lorentz-Maxwell  equations 
solved  self-consistently  in  weak  optical  fields.  The  integral  equation  is  formed  by  summing  over 
all  electrons  in  the  beam.  The  result  has  no  reference  to  the  electron  phases  and  only  involves 
the  optical  field  evolution.  High  gain  growth  rates  can  be  recovered  from  the  integral  equation 
analytically;  the  integral  can  be  solved  numerically  for  more  complicated  cases.  Contour  maps 
are  used  to  show  the  gain  degradation  due  to  an  electron  beam  energy  spread  and  an  electron 
beam  angular  spread.  In  the  limit  of  low  gain,  the  gain  spectrum  is  related  to  the  spontaneous 
emission  line-shape  through  successively  higher  derivatives.  In  the  limit  of  high  gain,  it  is  shown 
that  the  growth  rate  becomes  less  susceptible  to  degradation  from  the  electron  beam  quality.  The 
reference  is 

"Free-Electron  Laser  Gain  Degradation  and  Electron  Beam  Quality",  W.  B.  Colson,  J.  C. 

Gallardo  and  P.  M.  Bosco,  submitted  to  Physical  Review  A. 

The  document  report  is  included  at  the  end  of  this  report. 

Understanding  the  Free  Electron  Laser  Optical  Self-Guiding  Effect 

In  the  free  electron  laser  (FEL)  interaction  the  optical  light  wave  and  electrons  travel 
together  through  the  periodic  magnetic  field  of  the  undulator.  Without  the  electron  beam 
interaction,  the  optical  wavefronts  propagate  in  vacuum  and  diffract  according  to  Maxwell’s  wave 
equation.  Depending  on  the  wavelength,  diffraction  can  provide  an  important  limitation  to  the 
effective  interaction  length  of  the  undulator.  This  would  reduce  the  coupling  of  the  electrons  and 
light,  because  the  light  diffracts  away  the  smaller  electron  beam. 

An  important  FEL  effect  is  the  optical  guiding  of  light  which  focuses  the  light  wave  back  into 
the  electron  beam  so  the  interaction  strength  can  remain  strong  over  long  undulators.  The 
existence  of  this  effect  is  crucial  to  XUV  application  and  high  power  FELs.  This  brief  section 
contains  a  simple  understanding  of  how  the  self-guiding  effect  can  work  favorably  for  FELs. 


_# _ statement. 

1  do8x=l,n-i 

2  do8y=1,n-1 

3  artmp=ar(x,y) 

4  ardif=ar(x+l  ,y)+ar(x-1  ,y)+ar(x,y+1)+ar(x,y-1) 

5  aidif=ai(x+1  ,y)+ai(x-1  ,y)+ai(x,y+1  )+ai(x,y-1 ) 

6  ar(x,y)=ar(x,y)-C*(aidif-4*ai(x,y)) 

7  ai(x,y)=ai(x,y)+C*(ardif-4*artmp) 

8  continue 


The  first  two  steps  set  up  Fortran  loops  to  line  8  over  the  array  surface  Oct  <n  and  0<y  <n . 
The  arrays  a  =  ( ar.ai )  are  dimensioned  from  0  to  n.  In  order  to  avoid  sampling  array  elements 
outside  the  window  of  width  w,  the  loops  only  cover  sites  inside  the  edge.  The  third  step 
temporarily  saves  the  real  part  of  the  field  at  the  point  {x,y)  to  update  the  imaginary  part  in  step  7. 
Steps  4  and  5  evaluate  the  Laplacian  for  the  real  and  imaginary  parts  of  a(x,y)  using  the  nearest 
neighbor  differences  at  x  ±  1  and  y  ±  1.  The  real  part  of  a(x,y)  is  updated  in  step  6  with  the 
coefficient  C  =  Ax/4dx2  where  dx  =w/n.  The  number  of  time  steps  1/At,  and  the  window  width,  w , 
must  both  be  large  enough  so  that  C  is  a  small  number.  The  sampling  in  n  cannot  be  too  large. 
The  reason  for  this  is  that  diffraction  in  one  time  step  must  not  take  significant  light  past  nearest 
neighbors. 

We  see  that  the  important  small  quantity  in  this  expansion  is  Ax/4dx2.  This  algorithm  has 
been  used  on  an  IBM  PC  to  explore  diffraction  properties  in  the  free  electron  laser.  In  order  to 
understand  the  step  by  step  process  of  diffraction,  now  rewrite  the  diffraction  equations  for  one 
small  step  Ax  in  a  slightly  different  form. 

The  sum  of  surrounding  nearest  neighbor  points  divided  by  4  is  defined  as  the  average 
value  of  the  real  or  imaginary  parts  of  the  field  a.  The  step  equation  for  any  field  point  a{x,y)  is 
then  affected  by  its  nearest  neighbors  Aa=-i(a  -  a)  Axldx2.  The  optical  phase  is  given  by 
A<(>  =  -  (ar  -  ar)  Ax  /  \  a  \  dx2.  It  is  the  optical  phase  over  the  wavefront  that  determines  that 
diffraction  focuses  the  wave  into  or  out  of  the  the  electron  beam. 

A  wavefront  that  has  curvature  such  that  ar  <ar,  then  the  optical  phase  shift  will  be  positive. 
This  corresponds  to  focusing  of  the  wave.  If  ar  >  ar  then  the  phase  shift  will  be  negative  and  the 
wave  expands  in  size  away  from  the  electron  beam  reducing  coupling. 

The  FEL  interaction  also  changes  the  optical  phase  with  A$fel,  but  opposite  in  sign  to 
natural  diffraction.  Diffraction  causes  a  negative  phase  change  while  the  FEL  interaction  causes  a 
positive  phase  change.  The  FEL  interaction  therefore  decreases  the  effect  of  natural  diffraction 
and  focuses  the  wave  back  into  the  electron  beam. 
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1.  INTRODUCTION 

The  free  electron  laser  (FEL)  uses  a  high  quality  relativistic  beam  of  electrons  passing 
through  a  periodic  magnetic  field  to  amplify  a  co-propagating  optical  wave  (1-4).  In  an  oscillator 
configuration,  the  light  is  stored  between  the  mirrors  of  an  open  optical  resonator,  as  shown  in 
Figure  1.  In  an  amplifier  configuration,  the  optical  wave  and  an  intense  electron  beam  pass 
through  a  transversely  undulating  magnetic  field  to  achieve  high  gain.  In  either  case,  the 
electrons  must  spatially  overlap  the  optical  mode  for  good  coupling.  Typically,  the  peak  electron 
beam  current  varies  from  several  amperes  to  many  hundreds  of  amperes,  and  the  electron  energy 
ranges  from  a  few  MeV  to  a  few  GeV.  The  electrons  are  the  power  source  in  a  FEL,  and  provide 
from  a  megawatt  to  more  than  a  gigawatt  flowing  through  the  resonator  or  amplifier  system.  The 
undulator  resonantly  couples  the  electrons  to  the  transverse  electrical  field  of  the  optical  wave  in 


vacuum. 
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Figure  1.  The  basic  elements  of  a  free  electron  laser  (FEL)  oscillator  are  a  high  quality 
relativistic  electron  beam,  an  undulator  magnet  that  causes  the  electrons  to  wiggle,  and  the 
resonant  optical  cavity  to  provide  feedback. 

The  basic  mechanism  of  the  coherent  energy  exchange  is  the  bunching  of  the  electrons  at 
optical  wavelengths.  Since  the  power  source  is  large,  even  small  coupling  can  result  in  a  powerful 
laser.  Energy  extraction  of  5%  of  the  electron  beam  energy  has  already  been  demonstrated.  The 
electron  beam  quality  is  crucial  in  maintaining  the  coupling  over  a  significant  interaction  distance, 
and  of  central  importance  to  all  FEL  systems  is  the  magnetic  undulator.  The  peak  undulator  field 
strength  is  usually  several  kG  and  can  be  constructed  from  coil  windings  or  permanent  magnets. 
In  the  top  part  of  Figure  2,  the  Halbach  undulator  design  is  shown  for  one  period.  The  field  can  be 
achieve  1,  to  a  good  approximation,  using  permanent  magnets  made  out  of  rare  earth  compounds, 
a  technique  developed  by  K.  Halbach  (5)  and  now  employed  in  most  undulators.  The  undulator 
wavelength  is  in  the  range  of  a  few  centimeters  and  the  undulator  length  extends  for  a  few  meters, 
so  that  there  are  several  hundreds  periods  for  the  interaction  (6-8).  The  polarization  of  the 
undulator  can  be  either  linear  or  circular  or  a  combination  (9).  The  optical  wave  has  the  same 
polarization  as  the  undulator  driving  it.  This  is  an  illustration  of  the  FEL’s  most  important 
attribute--the  flexibility  of  its  design  characteristics. 

The  transverse  undulations  of  electrons  with  energy  ymc  2  generate  spontaneous  emission  in 
a  forward  cone  of  angular  width  y-1.  When  the  undulator  fields  are  strong  enough  so  that  the 
amplitude  of  the  cone’s  oscillation  off-axis  is  comparable  to  the  cone's  width,  a  detector  on-axis  at 
infinity  will  begin  to  see  several  radiation  harmonics  (10).  If  the  angular  deviations  of  the  cone  are 
larger,  then  the  spectrum  becomes  broad-band  like  the  synchrotron  emission  from  a  bending 
magnet.  The  total  emission  energy  from  a  bending  magnet  and  a  FEL  undulator  are  similar,  but 
the  FEL  spectrum  is  confined  to  a  relatively  narrow  bandwidth  because  the  electron  motion  is 
periodic  and  the  radiation  cone  stays  on  the  undulator  axis.  The  FEL  gain-bandwidth  falls  within 
the  narrow  spontaneous  emission  spectrum  that  is  determined  by  the  number  of  undulator 
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periods.  The  laser  line-width  can  be  much  narrower  than  the  spontaneous  line-width  as  in  an 
atomic  laser;  the  narrow  line  and  long  coherence  length  are  established  by  mode  competition. 


uniform  electrons  bunched  electrons 

Figure  2.  A  practical  design  for  constructing  the  undulator  field  is  shown  at  the  top  where 
eight  permanent  magnets  are  used  to  form  one  undulator  period.  The  arrows  show  the 
directions  of  the  permanent  magnetic  field.  The  interaction  of  an  initially  azimuthally  uniform 
electron  beam  with  the  radiation  in  a  FEL  causes  the  electron  beam  to  bunch  in  an  optical 
wavelength.  It  is  this  bunching  that  c"  es  coherent  radiation. 

The  laser  frequencies  driven  by  the  FEL  mechanism  are  much  higher  than  the  oscillation 
frequency  of  the  electrons  in  the  undulator.  This  is  due  to  a  large  Lorentz  contraction  of  the 
undulator  wavelength  and  a  large  relativistic  Doppler  shift  of  the  emitted  radiation  in  the  forward 
direction.  The  relation  between  the  undulator  wavelength  Xu,  Ihe  optical  wavelength  X,  and  the 
electron  beam  energy  is  then  X  =  Xul2y2  and  the  mechanism  can  be  described  as  stimulated 
Compton  backscattering.  It  is  the  relativistic  factor  2y 2  that  allows  the  FEL  to  reach  short 
wavelengths.  Low  energy  beams  (5  MeV)  are  being  used  to  reach  wavelengths  longer  than 
atomic  lasers  (500  pm)  and  high  energy  beams  (1  GeV)  are  used  for  x-rays  (500& ),  as  shown  in 
Table  1  (11-27).  The  FEL  system  is  also  continuously  tunable  merely  by  changing  the  electron 
energy  of  the  electron  source.  Figure  3  shows  some  FEL  system  configurations,  which  are 
explained  more  fully  in  Section  2. 


•V-\  s*  «* 
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Table  1 .  Operation  of  free  electron  lasers 


Name  (Ref.) 


USCB  (20) 

LLNL  (18) 
Hughes  (21) 
Erevan  (24) 
Novosibirsk  (27) 


Year  of  first 
operation 


Wavelength 


Peak  power 


Stanford  (11) 

1976,1977 

10  pm,  3.4  pm 

130  kW 

Columbia  (25) 

1977 

1.5  mm 

8  MW 

NRL  (26) 

1977 

400  pm 

1  MW 

NRL/Columbia  (23) 

1978 

400  pm 

1  MW 

LANL  (13) 

1981,82 

10.6  pm 

10  MW 

NRL  (12) 

1981 

4.6-3.1  mm 

75  MW 

1983 

35  GHz 

17  MW 

Orsay  (15) 

1981,1983 

6500  A 

60  mW 

MSNW  (14) 

1982 

10.6  pm 

b 

Frascati  (22) 

1983 

5145  A 

C 

TRW  (16) 

1983 

1 .57  pm 

1.2  MW 

A  =  amplifier,  O  =  oscillator,  ASE  =  amplified  spontaneous  emission. 
b  Output  power  not  measured,  but  peak  loss  of  electron  energy  was  observed  to  be  9%. 
c  With  an  input  laser  power  of  6  W,  a  gain  of  3x1 0"4  was  measured. 
d  A  gain  of  1 .5%  was  measured. 

Figure  2  illustrates  the  basic  bunching  mechanism  used  to  obtain  coherent  radiation.  The 
electrons  leaving  the  accelerator  are  randomly  positioned  over  many  optical  wavelengths.  There 
are  typically  107  electrons,  or  more,  in  each  section  of  the  electron  beam  one  optical  wavelength 
long.  As  the  light  and  electrons  interact  at  the  beginning  of  the  undulator,  some  electrons  gain 
energy  and  some  lose  energy.  Those  that  gain  energy  move  a  little  faster  longitudinally  and  those 
that  lose  energy  move  a  little  slower;  this  creates  one  bunch  in  each  optical  wavelength. 


Some  FEL  Configurations 


Single  Pass 
|faoc  |  It  If 
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Figure  3.  Free  electron  lasers  can  be  in  a  variety  of  configurations,  as  depicted  here.  In 
fact,  three  of  the  five  types  have  already  operated. 

FELs  have  been  described  in  a  number  of  articles  in  the  general  scientific  press  (28-37).  In 
addition,  there  is  a  textbook  and  a  number  of  review  articles  on  the  subject  (38-40)  and  two 
Special  issues  Of  IEEE  Journal  of  Quantum  Electronics  contained  many  papers  on  FELs  (41,  42). 
Finally,  there  are  six  volumes  of  conference  proceedings  that  contain  hundreds  of  papers  and 
provide  a  good  introduction  to  the  FEL  literature  (43-48). 

2.  GENERALITIES 

2.1  History 

The  historical  development  of  FELs  can  be  traced  back  to  the  microwave  tubes,  backward 
wave  oscillators,  traveling-wave  tubes,  magnetrons,  and  klystrons  of  the  1940s,  shown  at  the  top 
of  Figure  4.  The  traveling-wave  tubes  were  similar  in  structure  to  the  FEL  in  that  they  used  mildly 
relativistic  electrons  traveling  through  periodically  undulating  electric  or  magnetic  fields  inside  a 
wave  guide.  The  radiation  wavelengths  produced  were  in  the  centimeter  range.  A  characteristic 
of  all  such  devices  was  the  closed  structure  used  to  store  the  radiation.  These  electron  tubes 
were  tunable  by  changing  the  electron  energy  and  using  higher  harmonics,  and  high  efficiencies 
were  common.  While  the  Motz  (49)  tubes  used  the  same  configuration  as  the  FEL,  the  operating 
mechanism  was  different.  A  tube  that  used  the  same  mechanism  as  that  in  a  FEL  was  invented 


by  Phillips  (50),  but  J.M.J.  Madey,  the  inventor  of  the  FEL,  was  unaware  of  the  Phillips  tube, 
although  he  did  know  of  Motz’s  work.  Shorter  wavelengths  could  not  be  reached  because 
electrons  did  not  oscillate  fast  enough  and  the  closed  resonator  could  not  be  made  small  enough. 

GENERRL  HISTORY 
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Figure  4.  FELs  grew  out  of  the  development  of  electron  tubes  and  atomic  lasers.  They 

retain  some  of  the  good  qualities  of  both. 

Atomic  lasers  were  invented  in  the  1960s  and  made  use  of  two  new  concepts  (51):  excited 
electrons  in  the  bound  states  of  atoms  or  molecules  oscillated  rapidly  to  produce  optical  radiation 
and  this  radiation  was  stored  in  open  optical  resonator. 

J.M.J.  Madey's  conception  of  the  FEL  (1)  came  from  a  mixture  of  the  attributes  of  microwave 
tubes  and  the  atomic  laser:  the  Motz  undulator  and  the  optical  resonator.  The  relativistic  Lorentz 
contraction  and  the  Doppler  shift  produced  high  frequencies  from  the  slower  oscillations  of  the 
electrons  traveling  near  the  speed  of  light.  The  FEL  is  tunable  just  as  the  early  electron  tubes 
were  but  it  works  at  short  wavelengths. 

Independently,  R.  Palmer,  P.  Csonka,  and  K.  Robinson  were  working  on  the  coherent 
emission  of  radiation  by  relativistic  electron  beams  (52). 


2.2  Basic  Concepts 


A  good  theoretical  approach  to  FELs  is  to  solve  the  relativistic  particle  dynamics  and  couple 
the  solutions  to  the  optical  wave  equation.  The  more  sophisticated  analytical  methods  employed 
in  the  analysis  of  plasmas  and  lasers  are  appropriate,  but  generally  not  needed.  The  first  classical 
theory  was  introduced  by  M.O.  Scully,  F.  Hopf,  et  al  [53]. 

The  initial  electron  density  has  no  structure  on  the  scale  of  the  FEL  optical  wavelength. 
Individual  electrons  are  only  influenced  by  the  radiation  field,  the  undulator  magnetic  field,  and 
possibly  the  Coulomb  fields  of  other  electrons  if  the  density  is  large  enough.  For  typical  undulator 
fields  and  wavelengths,  the  radiation  emitted  spontaneously  after  just  one  pass  is  sufficient  to 
define  a  classical  wave.  The  Lorentz  force  equations  for  an  electron  are 

dt  me 


(ib) 

(ic) 


where  Z?r  and  Br  are  the  optical  electric  and  magnetic  fields,  Z?„  is  the  undulator  field,  e  =  \e\  is 
the  electron  charge  magnitude,  c  is  the  speed  of  light,  m  is  the  electron  mass,  is  the  electron 
velocity,  and  ymc2  is  the  electron  energy.  Only  four  of  the  five  equations  (Equations  1)  are  needed 
to  completely  specify  the  problem.  The  undulator  axis  is  taken  along  the  z  axis  so  that  the 
transverse  optical  force  with  contributions  from  both  Z?r  and  is  proportional  to  |Z?r  |  (1  -  pj.  For 
relativistic  electrons  (1-pj)=1/2y2,  so  that  the  transverse  optical  force  is  small;  the  optical  electric 
and  magnetic  forces  combine  almost  to  cancel  when  y  »  1 . 

In  order  to  couple  energy  out  of  the  electron  beam,  the  time  average  of  ft/?,  must  be 
nonzero  during  the  interaction  time  in  the  undulator.  The  role  of  the  undulator  is  to  rotate  the 
transverse  electron  velocity  as  the  field  Sr  passes  over  it.  Note  that  in  Equation  (la)  the 
transverse  electron  motion  is  determined  primarily  by  the  undulator  magnet  since  the  transverse 
optical  force  is  small.  However,  a  randomly  distributed  electron  beam  will  have  <  >  =  0  with 

no  net  energy  transfer.  But,  an  energy  modulation  alters  the  electron  z  velocities  to  cause 
bunching  and  coherent  emission.  While  deflections  off  the  mode  axis  are  necessary  for  coupling, 
they  cannot  be  too  large,  since  the  optical  mode  has  a  limited  radial  extent. 


A  suitable  undulator  field  (6)  around  the  mode  axis  is 

Bx  -  -B  { [  1  +  k2  (3x2  +  y  2)  /  8  ]  COS(*uz )  -  [  k2xy  /  4  ]  sin(Jfc„z ) } 

By-B  {[1+k2(x2  +  3y2)/8]s\n{kuz)-[k2xy  /  4  ]  cosfcz ) }  (2) 

B,~-B  ^  +  k2  (x2  +  y2)/8][xsin(^z)  +  ycos(^z)]  , 

where  B  is  the  peak  field  strength  and  ku  =  2n/ku  is  the  undulator  wavelength.  The  electron 
beams  suitable  for  FELs  must  be  sufficiently  aligned  so  that  the  transverse  excursions  are  small 
compared  to  K-  The  average  magnetic  field  strength  increases  off-axis  so  that  the  electrons  are 
focused  toward  the  axis.  When  electrons  are  focused  back  toward  the  undulator  axis,  the 
transverse  oscillations  are  called  betatron  oscillations.  Typical  transverse  excursions  are  small 
enough  that  kux  and  kuy  are  negligible. 

With  a  small,  high  quality  beam,  the  undulator  field  sampled  by  electrons  is 
(B cos(kuz),  B sin(jtMz  ),0)  and  the  orbits  which  are  helical,  are 

[  ( -KH)  cos (kuz),  ( -K  iy)s\n(kuz),  p,  ]  (3) 

where  P2  =  1— (1  +K2)l2y2  and  K  =eBXu/2nmc2.  Typically  K  =  1  and  one  sees  that  the  transverse 
oscillations  are  small. 

The  optical  field  polarization  that  best  couples  to  the  above  trajectory  is  given  by  the  vector 
potential 

E(t) 

X{z,t)  = - -  {sin[*z-co/+4>(r)], cos[fcz -cor +<(>(* )], 0}  ,  (4) 

k 

where  E(t)  is  the  electric  field  magnitude,  \  =  2ncl<a  =  2nlk  is  the  optical  carrier  wavelength,  and 
4>(/)  is  the  optical  phase.  No  x  ory  dependence  is  included  in  X ,  for  now,  since  we  assume  the 
electrons  remain  well  inside  the  optical  mode  waist.  The  optical  electric  field  is  Xr  =  -c~'dX/dt. 
Inserting/?,  and  Equation  (3)  into  Equation  (1b)  we  have 

dy  eKE 

—  =  -  cos[(*„  +  *)z  -at  +$]  .  (5) 

dt 

A  particularly  useful  form  of  Equation  (5)  may  be  obtained  in  the  case  where  the  fractional 
energy  change  &y/y«1.  Define  the  electron  phase  Z[t)  =  {ku+k)z(t)-u>t,  then  eliminate  y[t) 
from  Equation  (5)  to  get 


(6) 


- 2=~j~=  1°  lC0S(C  +  <i>)  . 

dx2  dx 

where  \a\  =4nNeKLE/y2mc2  is  the  dimensionless  optical  field  strength,  x  =  ci/L  is  the 
dimensionless  time,  L  =NXU  is  the  undulator  length  so  that  0<x<1,  and  v  =  dXjdx  is  the  electron 
phase  velocity.  The  electron  dynamics  have  been  put  in  the  form  of  a  pendulum  equation  [54], 

The  evolution  of  each  electron  entering  the  FEL  undulator  follows  Equation  (6).  Individual 
electrons  are  identified  by  their  initial  conditions  C(0)  =  Co  and  v{0)  =  v0  =  L[ku  +*)pi(0)-*].  In 
weak  fields  |  a  |«ti,  and  when  |  a  j»jc,  the  fields  are  considered  strong  because  the  phases 
evolve  significantly  in  the  time  t<  1.  Experiments  are  usually  designed  so  that  the  spread  in 
electron  velocities  does  not  cause  a  spread  in  v  0  greater  than  n.  This  can  be  adjusted  by  keeping 
the  length  L  small  enough,  but  a  better  beam  quality  allows  a  greater  length  L  and  much  more 


The  optical  wave  is  governed  by  the  wave  equation  driven  by  the  current  T,: 


2  2->2  J)  > 


where  the  {x,y)  dependence  has  been  dropped  (see  Section  5).  The  transverse  electron  current 
is  the  sum  of  all  particle  currents 

.  (8) 
m 

where  r*m(t)  is  the  trajectory  of  the  mth  electron  and  ^i  =  (pj,p>,0).  Even  the  spontaneous 
emission  spectrum  in  a  FEL  has  a  long  coherence  length  so  that  the  field  E(t)  and  phase  (J>(r)  can 
be  taken  to  vary  slowly  over  an  optical  period,  of1.  Then,  the  terms  containing  second  derivatives 
in  Equation  (7)  are  negligible  compared  to  terms  with  single  derivatives  and 

7 "(£«*)  =  -  =  -  * cKL  2  —  S(3)[J?-  K,(t) }  . 

dx  c  dx  m  Y 

Then,  the  wave  equation  has  the  simple  form 


««  .  -,r 

—  =  -j<e^>  , 

dx 


where  a  =  \a  \e‘*t  the  dimensionless  current  density  is  j  =8N(neKL)2p/y\c2,  p  is  the  electron 
particle  density,  and  the  angular  brackets  represent  a  normalized  average  over  the  electrons.  If 


electrons  are  bunched  at  the  phase  k,  then  the  optical  amplitude  is  driven  with  strength  j  during 
the  time  0<x<1  and  there  is  gain.  If  the  phase  jt/2  is  overpopulated,  then  the  optical  phase  <J> 
grows  with  little  gain.  Usually,  it  is  a  combination  of  |a  |  and  that  is  driven  because  the  electron 
bunching  is  not  perfect. 

Figure  5  shows  the  phase  space  evolution  of  a  periodic  section  of  the  electron  beam  in  the 
(C,v)  coordinates.  The  separatrix  path  shown  is  given  by  v/  =  2|a  | [1+sin(^  +  <j>)];  the  peak-to- 
peak  height  is  4 1  a  | 1/2  and  the  horizontal  position  is  determined  by  <J>.  The  “fluid"  of  electrons 
starts  equally  populating  all  phases  and  at  the  phase  velocity  v0  =  2.6  for  maximum  gain.  As  the 
electron  fluid  evolves  in  Figure  5  it  becomes  darker,  becoming  black  at  x  =  1.  The  final  bunching 
is  near  the  phase  n  and  the  gain  and  optical  phase  shift  evolution  are  shown  at  the  right.  The 
initial  optical  field  is  weak  a(0)  =  ao=1,  and  the  final  gain  determined  numerically  is 
G  =  [  |  a  ( 1 )  | 2  -  a  2(0)]/a 2(0)  =  0. 1 35; . 
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Figure  5.  The  electron  phase  space  follows  sample  electrons  through  the  undulator.  The 
separatrix  is  shown  as  a  guide  to  the  phase  space  paths.  The  electron  fluid  grows  darker  as 
it  passes  through  the  undulator.  (The  same  representation  is  employed  in  Figures  6  and  8.) 
Bunching  at  the  phase  k  leads  to  gain,  but  also  affects  the  optical  phase. 

While  we  have  made  a  few  assumptions,  the  “pendulum"  and  wave  equations  (Equations  6 
and  9)  form  a  simple,  powerful  description  of  the  FEL  (54).  They  are  valid  for  both  weak  ( |  a  |  «n) 
and  strong  (|«|»7i)  optical  fields  in  either  high  (j»  1)  or  low  (/« 1)  gain  conditions.  It  is 
generally  important  that  both  the  optical  field  amplitude  and  phase  are  included  in  the  description. 
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When  the  optica!  tields  are  weak,  Equations  (6)  and  (9)  can  be  easily  linearized  in  a  (x): 

a  =iy<expt-i(Co  +  v0x)Ki  >:  Ci=  l«  |cos(Co+v0x  +  <{>)  ,  (10) 


where  (')  =  </(  )ldi,  C=  Co  +  v’o't  +  Ci.  and  Ci  is  C  expanded  to  lowest  order  in  |a|.  For  a  uniform 
|  beam  distribution,  the  average  of  any  quantity  /  is  given  by 

</>=/  d^{f)l2K  , 

f, 

r, 

f’  the  electron  coordinates  can  be  removed  from  Equation  (10)  and  the  optical  field  is  determined  by 

■  the  roots  to  the  cubic  equation 

•> 


t 


L 


b 


af-iv0af-  ~  j  =  0  ,  (11) 

with  the  field  of  the  form 

3 

a  =  a  o  exp (— iv  0x)  £  c,  exp(a,x)  , 

r= 1 

where  the  cr  are  determined  by  initial  conditions.  If  |  v0|  »n  so  that  the  FEL  is  far  off-resonance, 
the  driving  term  j  is  negligible  and  the  trivia!  uninteresting  solution  a=a0  is  obtained,  i.e.,  no  gain. 
If  the  current  density  j  is  large,  so  that  v0  is  negligible,  the  important  real  root  is  ar  =  [ji2)v3^3l2 
giving  exponential  growth.  The  complex  field  is  then  described  by 
a  (x)  =  (a0/3)exp[  (;72)1/3(V3+i)x/2)  ],  and  the  gain  is  exponential  after  an  initial  bunching  time. 

Figure  6  shows  the  phase  space  evolution  in  the  high  gain  case  where  ;'  =  100.  The 
electrons  are  started  at  vo  =  0  to  show  how  gain  is  achieved  on-resonance.  Bunching  occurs  at 
the  phase  nl2,  but  in  the  high  gain  case  a  significant  optical  phase  shift  changes  the  position  of  the 
separatrix  so  that,  relative  to  the  optical  wave,  bunching  is  at  phase  rc.  The  resulting  exponential 
growth  and  phase  evolution  are  shown  on  the  right.  The  exponential  gain  only  occurs  after 
bunching  is  established. 


H 
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Figure  6.  In  the  high  gain  case,  there  is  a  substantial  optical  phase  change  shift,  which  then 
shifts  the  separatrix.  The  height  of  the  separatrix  is  proportional  to  the  \a\v2  and  grows 
with  the  high  gain. 

In  the  low  gain  case,  both  v0  and  j  are  important  in  Equation  (11).  The  gain  is  no  longer 
exponential  and  all  three  roots  are  needed  to  find  the  final  gain  at  x  =  1 ,  which  is  given  by 


[2-2cosv0-v0sinv0I 

G{v0)  =  -j - - - 

Vo 


j  ^d_  sin(vo/2) ' 
2  dv0  (vq/2) 


(12) 


The  gain  is  antisymmetric  in  v0  and  peaks  at  G  =  0.135/  with  v0  =  2.6.  Figure  7  shows  the  plot  of 
G  (v 0)  above  the  accompanying  optical  phase  shift  (}>(v0)  =  ;[  2sinv0  -  v0(l+cosv0)  ]  /v$.  Note  that 
the  gain  spectrum  can  be  written  as  the  derivative  of  the  spontaneous  emission  spectrum 
( sin(vo/2)  /  (vq/2)  ] 2  This  remains  true  for  a  large  class  of  undulator  designs  and  is  known  as  the 
Madey  theorem  [55].  The  theorem  states  that  when  an  undulator  design  produces  a  spectrum 
s  (v 0)  the  gain  is  proportional  to  the  slope  of  the  spectrum  ds(v0)ldv0.  A  second  theorem  relates 
the  "second  moment  of  the  mean  electron  energy  loss  evaluated  to  first  order  in  the  optical  field 
strength,”  <  [5-/1*] 2  >,  to  the  "mean  energy  loss  evaluated  to  second  order  in  the  optical  field 
strength,"  <  &/2)>  : 
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Figure  7.  The  final  gain  and  phase  of  the  optical  wave  are  plotted  as  a  function  of  v0. 
Experimental  points  are  superimposed  to  show  agreement  between  small  amplitude  theory 
and  experiment  (Orsay)  (89). 

In  the  FEL  oscillator,  gain  over  many  passes  leads  to  strong  fields.  The  spontaneous  fields 
either  experience  exponential  growth  or  the  repeated  gain  of  Equation  12.  In  stronger  fields 
where  \a\>K,  the  gain  process  changes  and  begins  to  depend  on  |«  |.  Electron  phases  now 
evolve  too  far  in  phase  space  and  bunching  is  difficult  to  maintain.  Figure  8  shows  electrons  in  a 
strong  field  a0=  |a(0)  |  =8.  The  separatrix  is  now  large  and  electrons  are  trapped  in  the  closed 
orbit  region  of  phase  space.  Those  near  the  harmonic  circular  paths  oscillate  around  the  phase 
tc/2  at  a  frequency  |  a  |  v2\  these  oscillations  are  called  synchrotron  oscillations.  There  is  a 
decrease  in  gain,  i.e.,  saturation.  When  the  gain  is  reduced  to  equal  the  FEL  system  losses, 
steady-state  operation  is  established. 
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Figure  8.  Phase  space  evolution  in  the  strong  field  regime.  The  "synchrotron"  motion  of  the 
particles  has  led  to  saturation  and  energy  is  no  longer  transferred  from  the  electrons  to  the 
optical  wave.  Even  in  saturation  the  phases  of  the  optical  wave  evolves. 

A  method  used  to  extend  the  saturation  limit  of  FELs  was  proposed  by  Kroll,  Morion  & 
Rosenbluth  (56)  and  is  called  the  tapered  undulator.  As  electrons  lose  energy  to  the  optical  wave, 
the  undulator  properties  can  be  modified  to  accommodate  the  new  electron  energy.  As  y 
decreases  either  the  modulator  wavelength,  X„,  or  field  strength  D ,  can  be  decreased  to  maintain 
resonance.  A  simple  case  is  that  in  which  both  B  and  change  along  the  undulator  so  that  K  is 
constant.  When  such  a  taper  is  included,  the  pendulum  equation  acquires  an  accelerating  term, 

5  =  L2dku(z)ldz , 

C  =  8+  \a  |  cos(C+(S>).  (13) 

In  the  absence  of  the  field  |  a  | ,  electrons  appear  to  be  "accelerated”  to  higher  phase  velocities.  In 
strong  fields,  about  half  the  electron  phases  are  trapped  near  the  phase  n,  which  drives  the  optical 
amplitude  and  gain.  Figure  9  shows  the  final  position  of  electrons  in  phase  space  after  trapping 
has  occurred  in  strong  fields  a0~  40  and  with  tapering  such  that  5  =  6rc.  The  untrapped  electrons 
are  seen  at  the  top  of  the  phase  space  picture  spread  over  the  phase  axis  randomly.  The  gain  is 
higher  than  would  be  possible  at  this  field  strength  without  tapering.  The  tapered  undulator  is  a 
good  example  of  the  design  flexibility  of  FELs.  The  undulator  structure  [length,  polarization, 


wavelength  profile,  field  profile  B(z),  etc]  are  all  features  that  can  be  modified  to  enhance 
performance  for  a  particular  FEL  application. 
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Figure  9.  In  a  tapered  FEL  some  electrons  are  trapped  near  the  phase  that  drives  the 
optical  wave.  The  untrapped  electrons  are  distributed  over  many  phases  and  do  not  drive 
the  wave. 

An  example  proposed  by  Vinokurov  &  Shrinsky  (57)  is  the  klystron  FEL  (sometimes  called  a 
transverse  optical  klystron  FEL,  or  TOK)  where  the  undulator  is  split  into  two  sections  separated 
by  a  drift  or  dispersive  section.  The  purpose  is  to  achieve  higher  gain  for  a  given  interaction 
length  L.  The  dispersive  section  acts  like  the  bending  magnet  of  an  electron  energy  analyzer. 
Small  variations  in  the  electron  phase  velocity  v  caused  by  the  first  undulator  section  are 
translated  into  phase  changes  A £  =  Dv  at  the  end  of  the  dispersive  magnet  and  the  parameter  D 
measures  the  strength  of  the  dispersive  field.  The  theoretical  description  of  the  field  and  the 
electrons  uses  Equations  6  and  9  with  A l,  =  Dv  applied  to  each  electron  at  t  =  1/2.  This  results  in 
a  higher  degree  of  bunching,  and  therefore  greater  gain  than  given  by  Equation  1 2. 

When  the  undulator  is  designed  to  have  linear  polarization,  only  the  definitions  of  variables 
in  Equations  6  and  9  change  while  the  form  of  the  equations  remains  the  same.  The  modifications 
are  «  ->  -.7,(4)  ],  j  -» ;'[70(^)  -J^)  ] 2,  where  ^  =  /cr2/2(1  +  K2),  and  B  becomes  the  rms 

undulator  field  strength. 


2.3  Transverse  Effects 

The  one-dimensional  analysis,  which  we  have  employed  up  to  this  point,  leaves  out  all 
transverse  effects  except  the  simple  periodic  undulator  motion. 

First  we  discuss  electron  beam  transverse  effects.  A  helical  undulator  provides  focusing  of 
the  electrons  in  both  transverse  planes.  Sometimes  a  longitudinal,  solenoidal  field  is  employed  so 
as  to  give  even  more  focusing.  For  some  devices  the  cyclotron  resonance  in  this  field  coincides, 
or  almost  coincides,  with  the  FEL  resonance  and  makes  the  interpretation  of  these  experiments 
more  complicated  (12).  On  the  other  hand  this  juxtaposition  appears  to  enhance  the  gain,  but  is 
limited  to  long  wavelength  applications  because  of  the  upper  limit  on  attainable  solenoidal  field 
strengths. 

For  planar  undulators  there  is  only  "natural”  focusing  in  the  plane  perpendicular  to  the 
sinusoidal  motion  and  the  betatron  wave  number  is  k^y  =  eBl^2mcy  in  the  non-wiggle  plane, 
where  B  is  the  peak  field.  The  resonance  condition  is  maintained  as  a  particle  undergoes 
betatron  oscillations.  In  the  wiggle  plane,  generally  some  focusing  is  required  (50,  58,  59). 
Quadrupoles,  although  they  give  focusing,  seriously  degrade  FEL  performance.  A  planar 
undulator  field  is 

fi>=-J5cosh(/tuy)  cos (kuz)  j^+fisinh^y)  sin(£uz)  t, 

so  that  the  motion  is 


,  dx  B 
x  = —  =  —  1  + 
dz  y  . 


sin(£Bz)  , 


and  hence  increases  as  y  increases.  This  increase  with  y  just  balances  the  decrease  of 
y'&dyldz  when  y  increases  and  causes  p,  to  be  constant.  E.  T.  Scharlemann  (60)  has  shown 
how  shaping  the  undulator  pole  faces  with  a  slight  parabolic  curvature  provides  horizontal 
focusing  while  maintaining  p2  a  constant  of  the  motion.  The  curvature  causes  the  field  to  increase 
off-axis  and  provides  focusing  in  both  x  and  y .  If  the  pole  face  is  given  by  y(x)  =  T0(1  -k?l4), 
then  the  focusing  will  be  the  same  in  x  and  y  and  the  electron  beam  cross  section  will  be  round. 

It  is  necessary,  in  any  real  FEL,  to  avoid  resonances  between  the  various  frequencies  to 
which  the  particles  are  subject.  For  example,  one  must  avoid  a  resonance  between  betatron 
oscillations  and  integral  multiples  of  \u.  Also,  one  must  avoid  the  usual  coupling  resonances 
between  the  betatron  oscillations  in  x  and  y .  There  is  another  kind  of  resonance  that  must  also 


be  avoided:  a  synchro-betatron  resonance  between  the  "synchrotron  motion"  of  trapped  electrons 
and  transverse  betatron  motion  (61 , 62). 

We  turn  now  to  transverse  effects  of  the  electromagnetic  wave.  The  simplest  effect  is  the 
excitation  of  cavity  modes  in  an  oscillator.  Figure  10  shows  this  phenomenon  in  a  computer 
simulation  of  the  original  Stanford  experiment  where  the  electron  beam  has  been  moved  off-axis 
to  excite  a  combination  of  higher  order  modes. 


Figure  10.  Typically  FELs  are  made  to  produce  the  fundamental  mode  in  an  optical 
resonator  with  a  Gaussian  shape  in  x  and  y .  A  higher  order  mode  is  excited  here  by  moving 
the  electron  beam  off  of  the  resonator  axis.  The  theoretical  calculation  employed  the 
parameters  of  the  original  Stanford  FEL  (11). 

The  Rayleigh  range  is  a  measure  of  the  effect  of  diffraction.  For  a  light  beam  of  radius  w, 
the  Rayleigh  range  z,  =  jw2  /  X  is  the  propagation  distance  over  which  the  optical  wavefront 
doubles  its  area.  In  a  proper  FEL  design  one  wants  good  overlap  between  the  electron  beam  and 
the  light  beam  over  the  whole  interaction  length  so  that  z,  should  be  comparable  to  L .  However,  if 
the  FEL  has  sufficiently  high  gain  it  can  provide  "guiding"  to  the  light  and  keep  it  within  the  electron 
beam  for  many  Rayleigh  lengths,  as  in  an  optical  fiber  (63,  64).  This  is  seen,  dramatically,  in 


Figure  1 1 .  Intense  electron  beams  going  through  a  FEL  can  provide  optical  guiding  of  the 
radiation.  In  the  absence  of  guiding,  the  radiation  would  diffract  out  of  the  electron  beam 
tong  before  the  end  of  the  undulator. 

A  FEL  provides  an  effective  index  of  refraction,  n ,  by  changing  the  optical  phase  along  the 
interaction  length: 

Re(«)-1-“^  =  7^—  <sin(C  +  4>)>  , 
k  dz  \a\Lk 


lm(n)  s 


1  d\a\  _ 
k  |  a  |  dz  1 1 


<cos(C +  <!>)>  • 


For  an  optical  fiber,  guiding  occurs  if  Re(V'2)  +  lm(K2)/2  >  1,  where  the  (complex)  fiver  parameter, 
V,  is  given  by  V2- (n2- 1)£>2*2,  where  b  is  the  electron  beam  radius.  Thus  one  can  readily 
determine  when  guiding  takes  place,  provided  one  can  evaluate  the  averages  over  particles  of 


sin(C  +  <|>)  and  cos(£  +  <j>).  When  there  is  gain,  we  know  that  the  averages  of  sin  and  cos  are 
nonzero. 

In  the  exponential  growth  regime  one  can  evaluate  the  averages  analytically  (63-65).  One 
simply  augments  the  wave  equation,  Equation  7,  with  V*  and  then  approximates  this  transverse 
derivative  with 


*r 

The  result  is  that  Equation  1 1  becomes  for  v0  =  0, 

ar  +  i  &?(L/zr)  -  (/ I2)j  =  0  , 

where  the  length  of  the  undulator  is  L .  Thus  the  effect  of  diffraction  and  optical  guiding  are 
included  in  a  one-dimensional  theory.  Extension  to  a  warm  beam  and  to  vo*0  can  be  found  in  the 
quoted  literature  (63-65). 

2.4  Longitudinal  Effects 

The  simple  pendulum  and  wave  equations  (Equations  6  and  9)  are  valid  for  a  single 
complex  field  a  =  \a  le'*  with  only  a  single  frequency,  the  carrier  frequency  co.  A  realistic  FEL 
oscillator,  or  amplifier,  produces  a  spectrum  of  frequencies  surrounding  the  carrier  wave.  Usually, 
the  coherence  length  extends  over  several  optical  wavelengths  so  that  the  slowly  varying 
amplitude  and  phase  approximation  remain  valid.  To  generalize  the  optical  field  representation  to 
many  modes,  the  single  complex  field  a  (x)  becomes  a  ( k  ,x)  or  a  (z  ,x). 

Driving  the  carrier  phase  $  in  the  center  of  the  optical  wavefront  will  focus  the  light  along  the 
electron  beam  path.  Even  in  low  gain,  diffraction  couples  the  transverse  and  longitudinal  waves. 
The  phase  profile  4>(z)  in  a  low  gain  oscillator  is  determined  by  the  resonator  mirrors  and  their 
Rayleigh  length  z .  This  causes  a  shift  in  frequency  and  a  shift  in  the  gain  spectrum  in  an  oscillator 
(66). 

Often,  the  lack  of  distinct  electron  energy  levels  leads  to  questions  about  the  ultimate 
coherence  capabilities  of  FELs.  In  both  the  FEL  and  atomic  laser,  a  long  coherence  length  and 
narrow  frequency  spectrum  is  determined  by  mode  competition,  not  by  energy  levels.  In  the  low 
gain  case,  the  weak  field  gain  per  pass  in  each  mode  is  given  by  Equation  12.  The  number  of 
modes  within  the  gain  bandwidth  is  about  y2  (typically  y 2 »  1).  Figure  12  shows  the  evolution  of 
100  optical  wavelengths,  around  resonance.  The  spontaneous  emission  above  resonance 


experiences  gain  on  every  pass,  while  other  wavelengths  receive  less  gain  or  absorption.  The 
vertical  scale  follows  the  photon  number  q(A)  over  six  orders  of  magnitude  in  100  passes.  The 
spectrum  clearly  narrows  as  mode  competition  continues.  The  photon  number  evolves  as 
exp [G  (A )np],  where  np  is  the  pass  number  in  the  low  gain  oscillator  where  modes  are  uncoupled. 

Figure  12.  The  growth  of  coherence  in  the 
optical  wave  is  shown  by  following  100  modes 
from  spontaneous  emission.  The  photon  density 
at  the  wavelengths  near  peak  gain  grows  more 
rapidly  than  the  densities  at  surrounding 
wavelengths.  This  narrows  the  spectrum  after 
only  100  passes.  Evidently  the  laser  can 
become  narrow-band. 


Short  pulse  effects  (67)  in  FELs  can  also  be  described  by  generalizing  the  field  to  a(k).  An 
essential  concept  is  "slippage";  this  is  the  distance  that  light  travels  over  the  electron  beam  while 
the  electrons  travel  through  the  undulator.  It  is  given  by  L(1  -  pz) -NX  using  the  FEL  resonance 
condition.  The  ratio  of  the  slippage  distance  NX  to  the  electron  pulse  length  determines 
whether  or  not  short  effects  are  important.  If  NX«c,,  then  the  pulse  is  considered  long,  and 
each  part  of  the  pulse  experiences  gain  proportional  to  the  local  density.  If  NX»  a,,  then  the  FEL 
has  short  pulses  and  the  modal  structure  of  the  pulse  is  comparable  to  the  gain  bandwidth,  ~N~\ 

Since  electrons  bunch  when  they  reach  the  trailing  edge  of  the  optical  pulse,  the  optical 
pulse  receives  more  gain  on  its  trailing  edge  than  on  its  leading  edge  and  behaves  as  if  it  is 
traveling  slower  than  the  speed  of  light,  c.  This  effect  is  called  "lethargy"  (68)  and  must  be 
considered  in  the  oscillator  FEL,  where  the  resonator  mirror  spacing  and  the  electron  pulse 
repetition  time  must  be  synchronized  (69,  70).  The  range  of  mirror  positions  to  achieve 
synchronism  is  astonishingly  small:  only  a  4-jjm  range  was  observed  in  the  Stanford  experiment. 
The  amount  of  synchronism  within  the  working  range  is  important  in  determining  the  laser  line- 
width  and  power. 

Other  longitudinal  effects  involve  long  pulses  in  the  FEL.  One  is  the  "trapped  particle" 
instability  analyzed  by  Kroll  &  Rosenbluth  (71).  The  synchrotron  frequency  |  a  |v2  can  mix  with 
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the  carrier  wave  and  produce  sideband  gain  in  the  FEL.  Figure  13  shows  the  growth  of  sideband 
structure  in  |  a  (z)  |  and  <j>(z).  A  window  section  of  a  long  pulse  is  four  slippage  distances  long 
(-2  < zlN\ <  2).  The  field  |a(z)|  is  plotted  at  the  top  left,  with  bright  regions  indicating  an 
intense  field  and  dark  regions  indicating  a  low  field  region.  The  pass  number  is  plotted  along  the 
vertical  axis.  The  "trapped  particle"  instability  starts  a  modulation  in  the  field  magnitude  |a(z)| 
and  the  phase  4>(z)  with  a  period  equal  to  the  slippage  distance.  The  final  spectrum,  the  Fourier 
transform  of  a(z),  is  shown  with  its  sideband  on  the  bottom  right;  above  is  the  weak  field  gain 
spectrum  for  reference.  The  final  electron  energy  spectrum  is  shown  above  the  gain  spectrum. 
The  power  and  net  gain  evolution  are  plotted  on  the  upper  right  as  a  function  of  pass  number  n . 
The  trapped  particle  instability  is  expected  in  nearly  all  FELs  that  saturate  because  of  strong 
fields. 
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Figure  13.  When  the  electron  synchrotron  oscillations  mix  with  the  carrier  wave,  sidebands 
can  be  formed.  Over  many  passes  the  optical  wave  develops  a  modulation  whose  period 
matches  the  synchrotron  period.  The  optical  power  increases  with  the  addition  of  energy  of 
the  sideband. 
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In  a  linearly  polarized  undulator,  the  electron  z  motion  is  more  complex  than  in  the  helical 
case  because  there  is  a  periodic  oscillation  of  the  electron  z  velocity  even  when  injected  perfectly. 
The  oscillation  in  z,Az,  is  given  by  khz  =  -  ^s\n(2kuct)  where  4  =  K2/2(1  +  k2).  Since  typically 
K~  1,  the  oscillations  are  a  sizeable  fraction  of  the  optical  carrier  wavelength  and  lead  to 
spontaneous  emission  and  gain  in  higher  optical  harmonics  (72).  To  generalize  Equations  6  and 
9  for  a  harmonic  hk ,  make  the  replacements  £  -»  h  £,  v  -*  hv . 

a  -*ah[J £)  -  7 {a+i)/2(7z £)  ],  and)  -*  jh[J (*_i)/a(A^)  —  J (*+iy2(A^)  ] 2 

The  form  of  the  equations  stays  the  same,  only  the  couplings  are  modified.  Note  that  there  is  gain 
only  in  the  odd  harmonics  h  =  1 ,3,5,  •  ■  • .  If  the  undulator  field  is  large  enough  so  that  K> 2,  then 
the  coupling  to  higher  harmonics  is  very  strong.  Several  of  the  FEL  experiments  to  date  have 
observed  coherent  emission  into  higher  harmonics,  and  it  should  prove  to  be  a  useful  technique 
for  reaching  shorter  wavelengths  in  a  FEL. 

3.  Free  Electron  Laser  Systems 

FELs  can  be  made  in  a  variety  of  configurations  as  is  depicted  schematically  in  Figure  3.  In 
Section  2  we  describe  in  more  detail  a  particular  linac  oscillator,  a  linac  amplifier,  and  a  storage 
ring  oscillator  experiment.  FEL  systems  are  rapidly  evolving  and  in  the  future  can  be  expected  to 
be  quite  different  from  those  described  here. 

In  Table  f  we  presented  a  compendium  of  those  FELs  that  have  operated.  Many  more  FEL 
devices  are  under  construction  and,  as  one  can  see  from  the  dates  in  Table  1 ,  these  devices  are 
being  brought  into  operation  at  an  ever-increasing  rate.  In  Table  2  (73-87)  we  present  a 
representative  list  of  FEL  accelerators. 

Of  great  importance  to  FELs  are  electron  beams  of  high  quality.  Two  figures  of  merit  of 
quality,  for  a  given  current,  are  energy  spread  and  brightness.  The  brightness  is  defined  by 
J  =tc2//y25V.  and  becomes  a  measure  of  "goodness,"  where  /  is  the  current  enclosed  within  the 
transverse  4-volume  (8V  ^  bxbx'byby').  For  uniform  phase  space  density,  the  brightness  can  be 
approximated  by  J  =  2J/y2bxbx'  byby'.  The  quality  of  a  beam  depends  upon  the  parameters  of 
the  accelerator,  the  type  of  accelerator,  and,  of  course,  with  what  care  it  is  aligned,  etc.  In  Table  2 
we  present  brightness  and  energy  spread  for  a  number  of  accelerators.  As  one  can  deduce,  the 
expected  performance  of  FELs  far  exceeds  present  achievements. 

The  development  of  FELs  has  been  the  result  of  both  theoretical  advances,  which  we 
emphasized  in  this  article,  and  of  experimental  advances.  In  fact,  without  the  latter,  we  would  only 
have  an  empty  theoretical  structure.  The  experimentalists  who  have  been  instrumental  in  the 


development  of  FELs  are  many  in  number  and,  of  course,  are  cited  in  the  references,  but  special 
note  should  be  taken  of  the  work  of  C.  A.  Brau,  D.  Prosnitz,  D.  A.  G.  Deacon,  J.  Eckstein,  L. 
Elias,  E.  Shaw,  S.  Skrinski,  B.  Kincaid,  C.  Pellegrini,  J.  M.  Ortega,  M.  W.  Poole.  A.  Renieri,  P. 
Elleaume,  T.  Smith,  A.  Gover,  J.  A.  Edighoffer,  J.  M.  Slater,  D.  Dattoli,  V.  Granatstein,  B. 
Newman,  R.  Warren,  T.  Marshall,  J.  Walsh,  R.  Panted,  J.  Pasour,  and  Y.  Petroff. 
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3.1  The  Linac  Oscillator 


The  experiment  of  the  TRW  Group  (15)  serves  to  illustrate  the  linac  oscillator.  The 
superconducting  accelerator  at  Stanford  has  a  bunch  length  4.2  ps,  a  peak  current  of  0.5  -  2.5  A, 
and  at  66  MeV  an  energy  spread  of  0.03%  and  a  beam  emittance  of  I.StcxIO-3  cm  •  rad.  The 
optical  cavity  had  mirrors  12.68  m  apart  with  a  7.5-m  radius  of  curvature.  At  the  optical 
wavelength  of  1.57  pm,  the  reflectivity  was  99.84%.  The  unduiator  consisted  of  pairs  of  linear 
arrays  of  SmC05  permanent  magnets  with  wavelength  Xu  =  3.6  cm  and  a  peak  field  of  2.9  kG. 

The  experiment  was  designed  to  study  the  effect  of  tapering.  Furthermore  they  devised  an 
optical  klystron  so  the  multicomponent  unduiator  had  the  following  structure.  First,  there  was  a 
prebuncher  section  of  15  periods,  then  a  magnetic  dispersion  section  of  two  periods  and  a  total 
length  of  58.6  cm.  Then  90  periods  followed  that  could  be  tapered  and,  finally,  15  periods  of 
constant  unduiator.  The  tapered  part  was  varied  to  be  a  0,  1 ,  and  2%  taper  in  energy.  Beam 
diagnostics  consisted  of  14  insertable  fluorescent  screens  so  as  to  be  sure  the  beam  was  steered 
properly  and  the  mirrors  were  aligned  using  a  green  light  laser. 

With  a  1%  taper,  the  FEL  had  an  average  output  laser  power  of  4W  and  the  peak  power  was 
1.2  MW.  Since  the  mirror  transmission  was  0.13%  on  each  end  of  the  cavity,  the  intracavity 
optical  power  was  1 1  GW  cm'2.  The  repetition  rate  was  10  Hz  and  the  macro-pulse  length  5  ms 
with  the  micropulse  of  4  ps.  The  radiation  fundamental  was  at  1.57pm  and  the  laser  bandwidth 
was  1 .3%. 

Above  threshold  for  the  laser,  the  power  increased  by  a  factor  of  1010  over  that  of  the 
spontaneous  radiation!  The  FEL  took  305  passages  at  a  gain  of  7%  per  pass  to  get  to  1 0%  of  the 
saturated  level.  The  experimenters  also  observed  coherent  radiation  at  the  second  and  third 
harmonic  of  1 .6  pm. 

A  study  was  made  of  the  effect  of  tapering  the  unduiator.  For  an  untapered  case  the 
electron  transfer  of  energy,  efficiency,  should  be  (1/2N).  The  efficiency  was  measured  to  be  0.4%, 
which  compares  well  with  the  expected  value.  With  a  1%  taper  the  electrons  clearly  divided  into 
two  groups:  trapped  and  untrapped.  Most  of  the  electrons  (60%)  were  trapped  and  decelerated 
1-1.8%  while  the  untrapped  electrons  were  unchanged  in  energy.  Thus  the  beneficial  effect  of 
tapering  was  demonstrated. 
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3.2  The  Linac  Amplifier 

The  experiment  of  the  LBL  /  LLNL  is  representative  of  linac  amplifier  FELs  (18.88).  The  FEL 
was  run  as  a  single-pass  amplifier  in  the  microwave  range  at  34.6  GHz.  The  input  signal  was 
supplied  by  a  magnetron  of  peak  power  60  kW  and  a  pulse  length  of  500  ns. 

Use  was  made  of  the  LLNL  Experimental  Test  Accelerator  (ETA)  (73)  to  provide  a  6-kA,  3.3 
MeV  beam  with  emittance  of  0.23  n  cm  •  rad.  An  emittance  filter  was  used  to  reduce  the  beam 
current  to  approximately  500  A  with  a  normalized  edge  emittance  of  0.47  n  cm  •  rad.  The  highly 
chromatic  transport  of  the  ETA  beam-line  and  matching  quadruples  results  in  a  15-ns,  nearly 
monoenergetic,  beam  delivered  to  the  interaction  region. 

The  undulator  magnet  was  three  meters  long,  and  the  undulator  period  was  9.8  cm.  The 
longitudinal  variation  of  the  undulator  field  provided  strong  vertical  focusing.  Horizontally  focusing 
quadrupole  magnets,  surrounding  the  undulator,  provided  horizontal  focusing  while  only  slightly 
reducing  the  vertical  focusing  and  negligibly  affecting  the  FEL  resonance  condition. 

The  interaction  waveguide  was  a  rectangular,  oversized  waveguide  immersed  in  the 
undulator.  The  inside  dimensions  of  the  waveguide  were  9.83  cm  wide  by  2.91  cm  high.  The 
electric  field  was  horizontal  and  coupled  to  the  TE0)  waveguide  mode,  which  was  excited  by  the 
input  microwave  signal. 

The  signal  gain  in  the  amplified  spontaneous  emission  mode  (no  microwave  input  signal) 
was  measured  and  it  was  found  that  the  microwave  signal  grew  at  a  rate  of  13.4  dB  per  meter  for 
a  beam  current  of  450  A.  Extrapolating  this  growth  back  to  the  origin,  one  finds  that  the  effective 
input  noise  was  0.35  W. 

The  amplifier  gain  was  studied  both  as  a  function  of  undulator  magnetic  field  intensity  and  as 
a  function  of  undulator  length.  The  peak  output  power  of  80  MW  was  achieved  for  both  the  2  m 
and  3  m  long  undulator.  The  amplifier  went  into  saturation  at  2.2  m;  beyond  this  point,  the 
amplified  output  power  first  decreased  and  then  near  3  m  started  to  increase  again.  The  gain  as  a 
function  of  undulator  length  showed  an  exponential  gain  of  approximately  15.6  dB  m"1  up  to 
saturation.  This  was  in  close  agreement  with  the  small  signal  gain  measurement.  The  gain 
curves  for  the  1-m  and  2-m  undulators  are  relatively  symmetric  about  the  peak,  while  the  gain 


curve  for  the  3  m  long  wiggler  shows  a  marked  asymmetry  with  a  plateau  on  the  long  wavelength 
side  of  the  curve.  This  asymmetry  at  saturation  is  also  shown  in  the  numerical  simulations. 

Study  of  excitation  of  other  than  the  TE0i  mode  and  study  of  the  effect  of  varying  the 
undulator  parameters  (so  as  to  avoid  saturation  at  80  MW)  are  to  be  undertaken  in  the  near  future. 
What  has  been  shown,  so  far,  is  that  a  FEL  can  be  operated  in  the  high  gain  regime  (Gain  > 
2500). 


s' 


33  Storage  Rings 

The  first,  and  so  far  the  only,  operation  of  a  storage  ring  FEL  oscillator  was  achieved  by  the 
Orsay-Stanford  collaboration  using  the  Orsay  ring  ACO  (16,  89).  This  laser  operated  in  the  visible 
range,  at  6500°,  and  produced  75  pW  average  power  or  60  mW  output  peak  power.  The 
intracavity  peak  optical  power  was  2  kW. 

The  ACO  storage  ring  has  a  circumference  of  22  m  and  was  operated  between  160  and  166 
MeV.  Two  bunches  were  employed,  with  the  average  current  between  16  and  100  mA.  The  rms 
bunch  length  was  (in  time  units)  0.5  to  1  ns  and  the  energy  spread  (rms)  0.9x10  “3  to  1.3X10-3. 
Because  of  the  strong  radiation  damping,  the  transverse  size  (rms)  was  0.3  to  0.5  mm, 
corresponding  to  an  angular  spread  of  0.1  to  0.2  mrad. 

The  optical  cavity  was  5.5  m  long  so  the  round-trip  time  resonated  with  the  1 1  m  between 
electron  bunches.  The  mirror  radius  was  3  m,  the  Rayleigh  range  1  m.  Although  the  mirror 
transmission  was  only  3x1 0~5,  the  round-trip  cavity  loss  was  7x1 (T4,  primarily  because  of 
absorption  in  the  mirror  dielectric.  In  fact,  there  was  mirror  degradation  due  to  the  radiation 
harmonics  of  the  undulator,  which  forced  the  experimentalists  to  operate  ACO  at  a  reduced 
energy  (originally  they  had  expected  to  be  at  240  MeV)  and  to  operate  the  undulator  at  reduced 
magnetic  field  (K  =  1.1  to  1.2),  both  effects  tending,  of  course,  to  reduce  the  flux  at  higher 
harmonics. 

The  permanent  magnet  undulator  had  17  periods  with  a  period  of  7.8  cm,  and  a  total  length 
of  1.33  m.  It  was  operated  as  an  optical  klystron  in  order  to  increase  the  gain  per  pass.  This 
increased  the  gain  by  a  factor  of  2  to  7  so  as  to  reach  2x1 0"4  per  pass.  Lasing  with  such  low  gain 
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required  careful  alignment  of  the  electron  beam  onto  the  axis  of  the  optical  cavity,  high  quality 
mirrors,  as  well  as  precise  synchronism  between  the  light  pulse  reflections  and  the  electron  bunch 
revolution  frequency.  The  detuning  curve  gave  only  a  1.6  pm  full  width  at  half  maximum  near 
laser  threshold. 

The  laser  time  pulse  structure  was  a  series  of  pulses  and  showed  the  electron  rf  synchrotron 
frequency  (13  kHz),  and  the  27.2-MHz  bunch  frequency.  The  time  sequence  of  pulses  is 
understood  as  a  consequence  of  theoretical  study  (90).  In  frequency  space  the  laser  had  three 
lines  (near  6500  A  )  with  the  dominant  one  at  6475  A.  All  the  lines,  corresponding  to  maximum 
gains  in  the  klystron  FEL,  were  in  the  TEMqo  mode.  The  width  of  the  lines  was  2-4  A.  Tunability 
was  over  150  A  and  limited  by  mirror  reflectivity. 

The  storage  ring  FEL  is  the  only  configuration  mentioned  where  the  FEL  teeds  back  on  the 
electron  source.  On  each  pass  the  working  FEL  "heats"  the  electron  beam  by  introducing  an 
energy  spread.  Synchrotron  radiation  P,yn  from  the  bending  magnets  in  the  ring  damps  the 
excitations.  The  laser  power  at  saturation  is  determined  by  thermodynamic  equilibrium,  which 
results  in  weak  fields;  this  is  the  Renieri  limit  (91),  =  Psyn/2N .  The  efficiency  of  the  FEL  was 

only  2.4x1 0~5,  which  is  0.4  of  the  prediction  of  Renieri  for  this  case. 

3.4  Extensions 

We  have  seen  that  FELs  can  be  expected  to  be  efficient,  powerful,  reliable,  tunable  sources 
of  radiation  in  a  wide  range  of  wavelengths.  In  fact,  FELs  have  already  been  made  to  operate 
from  the  microwave  range  down  to  the  visible  range.  It  is  reasonable  to  expect  that  soon  we  shall 
have  FELs  readily  available,  for  many  different  applications,  from  microwave  wavelengths  to  soft 
x-ray  wavelengths.  When  augmented  with  atomic  and  molecular  lasers  and  conventional  radio 
tube  sources,  we  can  expect  to  have  coherent  radiation  sources  throughout  the  radiation 
spectrum  (currently,  one  can  see  one’s  way  to  300  A  ). 

Why  then  should  one  develop  even  more  devices?  Clearly,  because  they  can  be  design  for 
special  purposes,  have  special  properties,  be  less  expensive,  more  efficient,  etc.  The 
development  of  FELs  is  far  from  completed  and  really  only  starting;  a  number  of  extensions  of 


FELs  appear  to  be  possible.  Here,  we  mention  a  few  of  them  and  refer  the  interested  reader  to 
the  appropriate  literature. 

In  the  microwave  range  it  is  possible  to  apply  a  longitudinal  magnetic  field  of  sufficient 
strength  that  the  cyclotron  frequency  resonates  with  the  radiation  frequency.  Thus  one  can 
arrange  a  device  where  there  is  coincidence  between  the  FEL  resonance  and  the  cyclotron 
resonance  as  described  in  Section  2.3  (12,  92,  93). 

It  is  possible  to  replace  the  undulator  with  an  electromagnetic  field.  The  attainable  magnetic 
field  of  an  rf  wave  is  less  than  that  of  a  static  or  pulsed  magnetic  field,  but  the  wavelength  of  the 
"undulator''  can  be  made  less  than  that  of  a  conventional  undulator.  Thus,  one  can  get  to  short 
wavelengths  with  a  low  energy  electron  beam.  The  use  of  an  rf  wave  as  an  undulator  has  already 
been  demonstrated  (94)  and  demonstration  has  been  made  of  an  electromagnetic  wave  undulator 
FEL  by  an  NRL  group  (95).  This  group  had  the  electron  beam  produce  500  MW  of  12.5-GHz 
radiation  through  a  backward  wave  oscillator  mechanism,  and  then  used  this  radiation  as  an 
undulator  for  FEL  action.  In  this  manner  they  produced  200-GHz  radiation  with  peak  power,  not 
yet  optimized,  of  0.35  MW.  The  Santa  Barbara  group  (96)  plans  to  employ  the  same  idea,  but  use 
the  FEL  mechanism  to  generate  the  rf  field  of  an  "undulator"  in  a  "two-stage  FEL." 

We  have  concentrated  upon  so-called  "Compton  regime  FELs”  where  there  is  a  strong 
interaction  between  the  electrons  and  the  optical  wave,  but  where  the  interaction  between 
electrons  is  small.  In  the  opposite  case,  where  the  electrons  interact  strongly  through  Coulomb 
forces  so  that  a  density  fluctuations,  or  plasmon,  description  of  the  electron  beam  is  more 
appropriate,  the  FEL  is  said  to  be  in  the  "Raman  regime."  An  understanding  of  the  collective 
regime,  the  Raman  regime,  is  more  difficult  than  that  of  the  Compton  regime  but  offers  distinctive 
features.  Experiments  (12)  have  demonstrated  6%  conversion  efficiency,  and  large  power 
emission  (75  MW)  in  this  regime.  One  can  expect  more  development  of  these  devices  in  future 
years  (97). 

An  interesting  extension  of  the  FEL  is  to  operate  in  a  dielectric  media  (98,  99).  Gas  loading, 
for  this  is  the  manner  proposed  to  realize  the  dielectric  media,  changes  the  phase-matching 
condition  and  so  allows  a  wider  parameter  space  than  the  vacuum  FEL.  In  fact,  this  extension 
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can  be  nontrivial  and  would  appear  to  allow  operation,  for  example,  at  smaller  undulator  magnetic 
fields  than  in  the  conventional  FEL.  The  resonance  condition,  for  relativistic  electrons  is 
n-1+X/Xu  =  (1+  K2)!2y2,  for  a  medium  having  an  index  of  refraction  n .  Note  that  the  (n  - 1 ) 
term  can  easily  be  comparable  to  the  usual  (1/2X2)  FEL  term.  One  can  think  of  this  device  as 
being  a  suitable  combination  of  the  Cerenkov  effect  and  the  FEL  resonance. 

Another  interesting  extension  of  a  conventional  FEL  is  to  have  an  undulator  in  an 
isochronous  storage  ring  (100,  101)  in  which  particles  with  different  energies  take  exactly  the 
same  time  to  go  around  the  ring.  Thus  bunching  at  optical  wavelengths  is  preserved  around  the 
ring.  Most  rings  do  not  have  this  property  and  thus  the  electron  bunch  on  entering  the  undulator  is 
essentially  a  "new  bunch"  with  random  phases.  Rings  can  be  made  isochronous,  to  some  degree, 
so  that  the  bunching  of  a  FEL  can  be  preserved.  Clearly  this  is  advantageous,  and  it  can  be  done 
so  as  preserve  far-infrared  wavelength  bunching  as  has  been  shown  on  BESSY  (102).  A  FEL 
using  this  concept  has  not  yet  been  made;  it  is  doubtful  that  the  technique  can  be  extended  into 
the  visible,  but  for  the  infrared  it  could  make  a  very  interesting  device. 

Finally,  it  should  be  emphasized  that  "pushing"  FELs  to  shorter  and  shorter  wavelengths,  as 
has  been  spearheaded  by  J.  M.  J.  Made''  and  C.  Pellegrini,  may  require  no  "new  inventions,”  but 
nevertheless  be  difficult  and  a  signifies! ,  extension.  This  subject,  as  one  might  expect,  has 
received  considerable  attention  (85,  87,  103,  104).  Suffice  it  to  say  here  that  it  appears  possible 
to  construct  a  FEL  oscillator  down  to  about  500  A,  and  a  single-pass  FEL  growing  from  noise  to 
about  300  A.  Just  what  the  limits  are  remains  to  be  seen,  but  extending  the  Orsay  achievement 
by  an  order  of  magnitude  appears  to  be  possible. 

ACKNOWLEDGEMENTS 

We  are  grateiul  for  support  from  the  US  Air  Force  Office  of  Scientific  Research  Grant  No. 
AFOSR-84-0079  (W.B.C.)  and  from  the  Office  of  Energy  Research,  US  Department  of  Energy 
under  contract  DE-AC03-76SF  00098  (A.M.S.). 


-31  - 


Literature  Cited 

1.  Madey,  J.M.J.  J.  Appl.  Phys.  42:1906  (1971). 

2.  Madey,  J.M.J.  Stimulated  emission  of  radiation  in  periodically  deflected  electron  beam.  US 
Patent  3,822,410  (1974). 

3.  Elias,  L.R.,  Fairbank,  W.M.,  Madey,  J.M.J.,  Schwettman,  H.A.,  Smith,  T.l.  Phys.  Rev.  Lett. 
36:717(1976). 

4.  Deacon,  D.A.G.,  et  al.  Phys.  Rev.  Lett.  38:892  (1977). 

5.  Halbach,  K.  See  Ref.  47,  p.  Cl -21 1. 

6.  Blewett,  J.P.,  Chasman,  R.  J.  Appl.  Phys.  48:2692  (1977). 

7.  Elias,  L.R.,  Madey,  J.M.J.  Rev.  Sci.  Instrum.  50:1339  (1979). 

8.  Diament,  P„  Phys.  Rev.  Lett.  A  23:2537  (1981). 

9.  Kim,  K.-J.  Nucl.  Instrum.  Method  219:425  (1984). 

10.  Colson,  W.B.  Free  electron  laser  theory.  PhD  dissertation,  Stanford  Univ.  (1977). 

1 1 .  Deacon,  D.A.G.,  et  al.  Phys.  Rev.  Lett.  38:892  (1977). 

12.  Gold,  S.M.,  et  al.  See  Ref.  48.  p.  350. 

13.  Newman,  B.E.,  et  al.  Phys.  Rev.  Lett.  38:118  (1977);  Newman,  B.E.,  Hohla,  K.,  Warren,  R., 
Goldstein,  J.C .IEEE  J.  Quantum  Electron.  QE17:1480  (1981 ). 

14.  Grossman,  W.M.,  et  al.  Phys.  Rev.  Lett.  38:52  (1977). 

15.  Edighoffer,  J.A.,  et  al.  Phys.  Rev.  Lett.  52:344  (1984). 

16.  Billardon,  M.,  et  al.  Phys.  Rev.  Lett.  51 :1652  (1983). 

17.  Pasour,  J.A.,  Lucey,  R.F.,  Kapetanakos,  C.A.  Phys.  Rev.  Lett.  53:1728  (1984). 

18.  Orzechowski,  T.  J.,  et  al.  Phys.  Rev.  Lett.  54:889  (1985). 

19.  Fajans,  J..  Bekefi,  G.,  Yin.  Y.Z.,  Lax,  B.  Phys.  Rev.  Lett.  53:  246  (1984). 

20.  Robinson,  A.L.  Reported  in  Science  226:154  (1984). 


• .  <■. 


• ,« . • .  * 
“>V‘ 

-y:« 

V\'V- 


SS38B 


•Srf 
«  •/«  ^ 


21.  Dolezal,  F„  Harvey,  R.,  Palmer,  A.  Paper  presented  at  9th  Ini.  Conf.  on  Infrared  and 
Millimeter  Waves,  Osaka  (Oct.  1984). 

22.  Barbini,  R.,  et  al.  See  Ref.  47,  p.  Cl -1. 

23.  McDermott,  B.D.,  et  al  Phys.  Rev.  Lett.  41:1368  (1978). 

24.  Omatuni,  T.S.,  Petrosian,  M.L.,  Petrosian,  B.V.,  Shahbazian,  N.T.,  Oganesian,  K.B.  Erevan 
EPI-727  (42) -84  (unpublished  1984). 

25.  Marshall,  T„  et  al.  Appl.  Phys.  Lett.  31 :320  (1977). 


nr."4 


a  '.y,-' 
/C'Xv 


-32- 

26.  Granatstein,  V.,  et  at.  Appl.  Phys.  Lett.  30:384  (1 977). 

27.  Artamonov,  et  al.  Nucl.  Instrum.  Methods  177:247  (1980);  Kornyuikhin,  G.A.  et  al.  Inst.  Nucl. 
Phys.,  Novosibirsk,  Internal  Report  (unpublihsed). 

28.  Rothenberg,  M.S.  Phys.  Today  29:17  (Feb.  1976). 

29.  ‘The  Free  Electron  Laser."  Sci.  Am.  236:62  (Jun.  1 977). 

30.  Schwarzschild,  B.  Phys.  Today  36:17  (Dec.  1 983). 

31 .  Lasers  and  Applications  2:40  (Oct.  1 983) . 

32.  Robinson,  A.  Science  221 :937  (Sep.  1983). 

33.  Robinson,  A.  Science  226:153  (Oct.  1984). 

34.  Schwarzschild,  B.  Phys.  Today  37:19  (Nov.  1984). 

35.  Schwarzschild,  B.  Phys.  Today  37:21  (Nov.  1984). 

36.  Robinson,  A.  Science  226:1300  (Dec.  1 984). 

37.  Thomsen,  D.E.  Sci.  News  1 26:359  (Dec.  1 984). 

38.  Marshall,  T.C.  Free  Electron  Lasers.  New  York:  MacMillan  (1 985). 

39.  Prosnitz,  D.  "Free  electron  lasers."  In  CRC  Handbook  of  Laser  Science  and  Technology  1,  ed. 
M.J.  Weber,  p.  425.  Boca  Raton:  CRC  Press  (1982). 

40.  Pellegrini,  C.  IEEE  Trans.  Nucl.  Sci.  NS-26:3791  (1979);  Morton,  P.L.  IEEE  Trans.  Nucl.  Sci. 
NS-28:3125  (1981). 

41 .  IEEE  J.  Quantum  Electron..  Vol.  QE-1 7,  No.  8  (Aug.  1981). 

42.  IEEE  J.  Quantum  Electron..  Vol.  QE-1 9.  No.  3  (Mar.  1 983). 

43.  Jacobs,  S.F.,  Sargent,  M.  Ill,  Scully,  M.O.,  eds.  Novel  sources  of  coherent  radiation.  In 
Physics  of  Quantum  Electronics,  Vol.  5,  Reading,  Mass:  Addison-Wesley  (1978). 

44.  Jacobs,  S.F.,  Pilloff,  H.S.,  Sargent,  M.  Ill,  Scully,  M.O.,  Spitzer,  R.,  eds.  Free  electron 
generators  of  coherent  radiation.  In  Physics  of  Quantum  Electronics,  Vol.  7,  Reading,  Mass: 
Addison-Wesley  (1980). 

45.  Jacobs,  S.F.,  Pilloff,  H.S.,  Sargent,  M.  Ill,  Scully,  M.O.,  Spitzer,  R.,  eds.  Free  electron 
generators  of  coherent  radiation.  In  Physics  of  Quantutn  Electronics,  Vol.  8,  Reading,  Mass: 
Addison-Wesley  (1982). 

46.  Jacobs,  S.F.,  Pilloff,  H.S.,  Sargent,  M.  Ill,  Scully,  M.O.,  Spitzer,  R.,  eds.  Free  electron 
generators  of  coherent  radiation.  In  Physics  of  Quantum  Electronics,  Vol.  9,  Reading,  Mass: 
Addison-Wesley  (1982). 

47.  Bendor  Free  Electron  Laser  Conference,  J.  Phys.  Colloq.  Cl ,  Suppl.  2,  p.  44  (Feb.  1983). 

48.  Brau,  C.A.,  Jacobs,  S.F.,  Scully,  M.O.,  eds.  Free  electron  generators  of  coherent  radiation, 
SPIE  (Int.  Soc.  Opt.  Eng.)  453  (1984). 


& 

_ 


.«  v  v 
-v  ;v.'v 

mi  ; 


49.  Motz,  H.  Thon,  W..  Whitehurst,  R.N.y.  Appl.  Phys.  24:826  (1953). 

50.  Phillips,  R.  M.  IRE  Trans.  Electron  Devices  1 7:231  (1960). 

51.  Schawlow,  A.L.,  Townes,  C.H.  Phys.  Rev.  112:1940  (1958);  Maiman,  T.H.  Phys.  Rev.  Lett. 
4:564  (1960). 

52.  Palmer,  R.B.  J.  Appl.  Phys.  43:3014  (1972);  Csonka,  P.  Part.  Accel.  8:225  (1978);  Robinson, 
K.,  unpublished  notes. 

53.  Hopf,  F.A.,  Meystre,  P.,  Scully,  M.O.,  Louisell,  W.H.  Phys.  Rev.  Lett.  37:1342  (1976). 

54.  Colson,  W.B.  Phys.  Lett.  64A:190  (1977);  Colson,  W.B.,  Ride,  S.K.  Phys.  Lett.  76A:379 
(1980);  Tang,  C.-M.,  Sprangle,  P.  Appl.  Phys.  53:831  (1982);  Shih,  C.-C.,  Yariv,  A.  IEEE  J. 
Quantum  Electron.  QE-1 7:1387  (1981);  Dattoli,  G.,  et  al.  IEEE  J.  Quantum  Electron.  QE- 
17:1371  (1981);  Bonifacio,  R.,  Pellegrini,  C.,  Narducci,  L.U.Opt.  Common.  50:373  (1984). 

55.  Madey,  J.M.J.  Nuovo  CimentoB  50:64  (1979). 

56.  Kroll,  N.M.,  Morton,  P.L.,  Rosenbluth,  M.N.  IEEE  J.  Quantum  Electron.  QE-27:1436-68 
(1981). 

57.  Vinokurov,  N.A.,  Shrinsky,  A.N.  Inst.  Nucl.  Phys.,  Novosibirsk,  Rep.  No.  INP  77-59 
(unpublished)  (1977). 

58.  Luccio,  A.,  Krinsky,  S.  See  Ref.  45,  p.  181  (1982). 

59.  Dattoli,  G.,  Renieri,  A.  Experimental  and  theoretical  aspects  of  the  free  electron  laser.  In 
Laser  Handbook  4,  ed.  M.L.  Stitch,  M.S.  Bass.  Amsterdam:  North  Holland  (To  be  published). 

60.  Scharlemann,  E.T.  Lawrence  Livermore  National  Laboratory  ELF  Note  105P  (published)  (1984). 

61 .  Kroll,  N.,  Rosenbluth,  M.  See  Ref.  44,  p.  147. 

62.  Fawley,  W.M.,  Prosnitz,  D.,  Scharlemann,  E.T.  submitted  Phys.  Rev.  A  30:1472  (1984); 
Lawrence  Livermore  Natl.  Lab.  Rep.  UCRL90838  (unpublished)  (1984). 

63.  Scharlemann,  E.T.,  Sessler,  A.M.,  Wurtele,  J.S.  Proc.  Como  Conf.,  Nucl.  Instrum.  Methods.  In 
press;  Lawrence  Livermore  Natl.  Lab.  Rep.  UCRL-9I476  (unpublished)  (1984). 

64.  Scharlemann,  E.T.,  Sessler,  A.M.,  Wurtele,  J.S.  Phys.  Rev.  Lett.  54:1925  (1985). 

65.  Moore,  G.T.  Opt.  Commun.  52:46  (1984);  Moore,  G.T.  Proc.  Como.  Conf.,  Nucl.  Instrum. 
Methods.  In  press;  IMO  Rep.  84-14  (unpublished)  (1984). 

66.  Colson,  W.B.,  Elleaume,  P.  Appl.  Phys.  B  29:1-9  (1982). 

67.  Colson,  W.B.,  Renieri,  A.  See  Ref.  47,  pp.  C1-11. 

68.  Al-Abawi,  H.,  Hopf,  F.A.,  Moore,  G.T.,  Scully,  M.O.  Opt.  Commun.  30:235  (1979). 

69.  Eckstein,  J.N.,  et  al.  See  Ref.  45,  p.  49. 

70.  Benson,  S.,  Deacon,  D.A.G.  et  al.  Phys.  Rev.  Lett.  48:235  (1982). 

71.  Kroll,  N.M.,  Rosenbluth,  M.N.  See  Ref.  44,  p.  147. 


-34- 


72.  Colson,  W.B.  IEEE  J.  Quantum  Electron.  QE-17:1417-27  (1981). 

73.  Fessenden,  T.J.,  et  al.  In  Proc.  4th  Int.  Conf.  on  High  Power  Electron  and  Ion  Beam  Research 
and  Technology,  ed.  H.J.  Doucet,  J.M.  Buzzi,  p.  813,  Palaiseau,  France:Ecole  Polytech. 

74.  Briggs,  R.J.,  et  al.  IEEE  Trans.  Nucl.  Sci.  NS-28.3360  (1981). 

75.  Takeda,  S.  Workshop  on  Generation  of  High  Fields,  Frascati,  Sep.  25-  Oct.  1,  1984.  Nucl. 
Instrum.  Methods.  In  press  (1985). 

76.  Warren,  R.M.,  et  al.  Lasers  84  Tech.  Dig.,  p.  20  (1984). 

77.  Elias,  L.  Lasers  84  Tech.  Dig.,  p.  25  (1984). 

78.  Boehmer,  H.,  et  al.  Phys.  Rev.  Lett.  48:141  (1982). 

79.  Birx,  D.L.,  et  al.  A  multipurpose  5-MeV  linear  induction  accelerator.  In  Proc.  IEEE  Power 
Modulator  Symp.,  June  1984;  Lawrence  Livermore  Natl.  Lab.  Rep.  UCRL-90554  (unpublished) 
(1984). 

80.  Shaw,  E.D.,  Chichester,  R.J.,  Sprenger,  W.C.  PLasers  84,  Tech.  Dig.,  p.  19  (1984). 

31 .  Smith,  S.D.,  et  al.  See  Ref.  45,  p.  275. 

82.  Bizzarri,  U.,  et  al.  Lasers 84  Tech.  Dig.,  p.  21  (1984). 

83.  Leiss.  J.,  Norris,  N.J.,  Wilson,  M.  A.  Part.  Accel.  10:223  (1980). 

84.  Billardon,  M.,  et  al.  See  Ref.  15,  p.  1652. 

85.  Deacon,  D.A.G.  Lasers  84  Tech.  Dig.,  p.  27  (1984). 

86.  Ballili,  J.P.,  et  al.  Report  of  Prospectives  for  Super  AGO,  LURE,  Orsay,  Internal  Report 
(unpublished)  (April  1981). 

87.  Peterson,  J.M.,  et  al.  Contributed  paper  to  the  Castelgandolfo  Conf.  Instrum.  Methods  in 
Physics  Research.  To  be  published  (1 985). 

88.  Orzechowski.  T.J.,  et  al.  See  Ref.  48,  p.  65. 

89.  Billardon,  M.,  et  al.  See  Ref.  47,  pp.  Cl  -29. 

90.  Elleaume,  P.  J.  Phys.  45:997  (1984). 

91 .  Renieri,  A.  Nuovo  Cimento  B  59:64  (1979). 

92.  Freund.  H.P.  See  Ref.  48,  p.  361;  Bernstein,  I.B.,  Friedland,  L.  Phys.  Rev.  A  23:816(1981). 

93.  Lin,  A.T.  See  Ref.  48,  p.  336;  Friedland,  L.,  Hirschfield,  J.L.  Phys.  Rev.  Lett.  44:1456  (1980). 

94.  Tsumoru,  S.,  et  al.  Jpn.  J.  Appl.  Phys.  22:844  (1983). 

95.  Granatstein,  V.L.,  Carmel,  Y.,  Gover,  A.  See  Ref.  48,  p.  344. 

96.  Elias,  L.B.  Phys.  Rev.  Lett.  42:977  (1980). 

97.  Freund,  H.P.,  Ganguly,  A.K.  See  Ref.  48,  p.  367. 

98.  Kimura,  W„  Wong,  D.,  Piestrup,  M„  Fauchet,  A.,  Edighoffer,  J.,  Pantell,  R.  IEEE  J.  Quantum 
Electron.  QE-18:239  (1982);  Walsh,  J.,  Johnson,  B.  See  Ref.  48,  p.  376. 


g 


-35- 


99.  Fauchet,  A.M.,  et  al.  See  Ref.  48,  p.  423. 

100.  Deacon,  D.A.G.  Phys.  Rev.  Lett.  44:449  (1980). 

101.  Deacon,  D.A.G.  Phys.  Rep.  76:349  (1981). 

102.  Gaupp,  A.  See  Ref.  47,  pp.  Cl -147. 

103.  Madey,  J.M.J.,  Pellegrini,  C.,  eds.  AIP  Conf.  Proc.,  Vol.  118,  Brookhaven,  Sep.  1983,  New 
York  (1984). 

104.  Proc.  Conf.  on  Coherent  and  Collective  Properties  in  the  Interaction  of  Relativistic  Electrons 
and  Electromagnetic  Radiation,  Como,  Sep.  1984.  Nucl.  Instrum.  Methods.  In  press  (1985). 


r»  [jy 

&& 

\  N  *, 

•v  - 

V/\ 

V-V-j 

V>. 

> 


SECURITY  CLASSIFICATION  OF  THIS  PAGE 


la.  REPORT  SECURITY  CLASSIFICATION 


REPORT  DOCUMENTATION  PAGE 


2a.  SECURITY  CLASSIFICATION  AUTHORITY 


2b.  DECLASSIFICATION /DOWNGRADING  SCHEDULE 


4.  PERFORMING  ORGANIZATION  REPORT  NUM8ER(S) 


1b.  RESTRICTIVE  MARKINGS 


3.  DISTRIBUTION/AVAILABILITY  OF  REPORT 


S.  MONITORING  ORGANIZATION  REPORT  NUMBER(S) 


6a.  NAME  OF  PERFORMING  ORGANIZATION 
Berkeley  Research  Associates 


6c.  ADDRESS  (City.  State,  and  ZIP  Code) 
P.0.  Box  241 
Berkeley,  CA  94701 


8a.  NAME  OF  FUNDING /SPONSORING 
ORGANIZATION 


8c.  ADDRESS  (City,  State,  and  ZIP  Code) 


6b.  OFFICE  SYMBOL  7a.  NAME  OF  MONITORING  ORGANIZATION 

(If  applicable)  Air  porce  office  of  Scientific  Research 


7b.  ADDRESS  (City,  State,  and  ZIP  Code) 
Bolling  Air  Force  Base 
Washington,  DC  20332 


8b.  OFFICE  SYMBOL  9.  PROCUREMENT  INSTRUMENT  IDENTIFICATION  NUMBER 
(If  applicable) 

F49620-85-C-0087 


10.  SOURCE  OF  FUNDING  NUMBERS 


PROGRAM 
ELEMENT  NO. 


PROJECT 

NO. 


WORK  UNIT 
ACCESSION  NO. 


1 1 .  TITLE  (Include  Security  Classification) 

The  Trapped-Particle  Instability  in  Free  Electron  Laser  Oscillators  and  Amplifiers 


12.  PERSONAL  AUTHOR(S) 

■JW.  B.  Colson 


13a.  TYPE  OF  REPORT 

Final  Scientific 


13b.  TIME  COVERED 
FROM  TO 


14.  DATE  OF  REPORT  (Year,  Month,  Day)  Il5  PAGE  COUNT 


17. _  COSATI  CODES _  18.  SUBJECT  TERMS  (Continue  on  reverse  if  necessary  and  identify  by  block  number) 

FIELD  I  group  |  SU8-GROUP  free  electron  lasers  laser  oscillator 


high-power  laser  laser  amplifier 

trapped-particle  instability 


•  19.  ABSTRACT  (Continue  on  reverse  if  necessary  and  identify  by  block  number) 

‘The  electrons  in  a  high-power  f ree-electron  laser  can  become  trapped  and  oscillate  in  deep 
potential  wells.  The  oscillating  component  of  the  current  drives  the  optical  field  at 
requencies  around  the  fundamental.  General  features  of  the  trapped-particle  instability 
f ree-eelctron  laser  oscillators  and  amplifiers  are  discussed.  Dimensionless  parameters 
lclarify  the  trends  for  a  wide  range  of  FEL  designs. 


I 


20  DISTRIBUTION /AVAILABILITY  OF  ABSTRACT 

□  UNCLASSIFIED/UNLIMITEO  □  SAME  AS  RPT.  □  DTIC  USERS 


21.  ABSTRACT  SECURITY  CLASSIFICATION 


22a.  NAME  OF  RESPONSIBLE  INDIVIDUAL 


DD  FORM  1473, 84  mar 


22b.  TELEPHONE  (Include  AreaCode)  I  22c.  OFFICE  SYMBOL 


83  APR  edition  may  be  used  until  exhausted. 
All  other  editions  are  obsolete. 


SECURITY  CLASSIFICATION  OF  THIS  PAGE 


y V. •y'y-\r\  ■  V  VV-  '■y^y''y\'-'y'~y-yy-y---:sy^'W.<\S'y:y~'y^y'y-y'.'-  VVv'-A 


The  Trapped-Particle  Instability  in 
Free  Electron  Laser  Oscillators  and  Amplifiers 

W.  B.  Colson 

Berkeley  Research  Associates,  P.O.  Box  241,  Berkeley,  California  94701 

ABSTRACT 

The  electrons  in  a  high-power  free-electron  laser  can  become  trapped  and 
oscillate  in  deep  potential  wells.  The  oscillating  component  of  the  current  drives 
the  optical  field  at  frequencies  around  the  fundamental.  General  features  of  the 
trapped-particle  instability  in  free-electron  laser  oscillators  and  amplifiers  are 
discussed.  Dimensionless  parameters  clarify  the  trends  for  a  wide  range  of  FEL 
designs. 

1.  Introduction 

The  free-electron  laser  (FEL)  uses  a  relativistic,  high-current  electron  beam  to  amplify  a 
copropagating  electromagnetic  wave  [1],  The  optical  wave  and  electrons  are  coupled  as  they 
pass  through  a  periodic,  transverse  magnetic  field  undulator.  A  common  goal  of  both  the  FEL 
oscillator  [2],  and  the  FEL  amplifier  [3]  is  large  optical  power.  In  strong  optical  fields,  the  electrons 
can  be  trapped  in  deep  potential  wells  formed  by  the  combined  optical  and  undulator  field  forces. 
Electrons  near  the  bottom  of  the  well  oscillate  in  harmonic  orbits  at  the  synchrotron  frequency. 
The  synchrotron  frequency  mixes  with  the  optical  carrier  wave  frequency  to  form  sidebands,  and 
the  growth  of  the  sideband  power  is  the  trapped-particle  instability.  In  both  the  amplifier  and  the 
overmoded  oscillator  configurations,  the  optical  amplitude  and  phase  are  free  to  evolve  over  a 
wide  range  so  that  complicated  optical  spectra  can  result  from  sideband  mixing.  In  this  paper,  we 
try  to  choose  examples  that  clearly  show  the  features  of  the  trapped-particle  instability  with  a 
minimum  of  complication  from  further  sideband  mixing. 

Historically,  the  instability  was  first  predicted  for  high-power,  low  gain  tapered  FEI.s  with 
many  synchrotron  oscillations  in  a  single  pass  through  the  undulator  [4],  Later,  it  was  seen  in  the 
simulations  of  pulsed  FEL  oscillators,  and  found  to  have  some  different  properties  than  originally 
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predicted  [5-7].  The  instability  has  been  clearly  observed  in  the  high-power  FEL  oscillator  at 
LANL  [2],  and  possibly  in  the  TRW/Stanford  oscillator  [8].  It  has  been  termed  the  sideband,  Kroll- 
Rosenbluth,  synchrotron,  Raman,  and  the  trapped-particte  instability  [9-23].  There  can  be  a  rich 
assortment  of  varied  effects  ranging  from  a  single  sideband  to  a  chaotic  frequency  spectrum  [22], 
The  shape  of  ultra-short  optical  pulse  can  be  controlled  to  form  shorter  subpulses  in  FEL 
oscillators  [5,6,13,21]. 

Any  FEL  working  at  high  power  is  susceptible  to  the  instability.  Tapered  undulators  are 
used  to  extend  the  FEL  performance  to  higher  power  levels  [9],  but  appear  to  have  less  sideband 
growth.  The  tapered  undulator  design  is  changed  along  its  length  to  maintain  resonance  with  the 
electrons  as  they  lose  energy  to  the  optical  wave.  We  review  simulations  of  the  instability  in  FEL 
oscillators,  and  present  new  results  of  the  instability  in  FEL  amplifiers.  Dimensionless  variables 
are  used  to  determine  important  regimes  of  operation. 

2.  Multimode  FEL  Theory 

The  wave  equation  determines  the  evolution  of  the  radiation  field  in  the  presence  of  the 
electron  current.  The  carrier  wave  with  frequency  co  =  /:c  =  2nc/\  has  the  single-mode  phase 
e,[a~m)  and  a  complex  slowly-varying  coefficient.  Transverse  (X,Y)  dependence  is  not  included, 
since  most  FEL  designs  attempt  to  minimize  the  effects  of  optical  diffraction  and  mode  coupling. 
The  intention  here  is  to  clarify  longitudinal  multimode  effects  without  distractions  from  other  FEL 
topics;  many  of  the  effects  can  be  approximately  characterized  by  adjusting  the  coupling  factors  in 
the  wave  equation.  Multimode  features  are  followed  with  a  complex  field  envelope 
a{z)  =  \a(z)  |e'w,)  evaluated  at  many  discrete  sites  z.  With  the  slowly-varying  amplitude  and 
phase  approximation  and  the  coordinate  change  Z->Z+ct,  the  wave  operator  reduces  to  a  single 
derivative  in  time  [24],  The  dimensionless  optical  field  envelope  is  a(z)  =  AnNcKLE{z)iylmcz 
where  N  is  the  number  of  undulator  periods,  e  is  the  electron  charge  magnitude,  L=NXq  =  2kNIIc0 
is  the  undulator  length,  K  =  eBXoi'Znmc2,  B  is  the  rms  undulator  field  strength,  m  is  the  electron 
mass,  c  is  the  speed  of  light,  E(z)  is  the  complex  optical  electric  field,  and  y 0mc2  is  the  resonant 
electron  energy,  =  k  (1+^2)/2L0  for  y0»1 .  A  resonant  electron  passes  through  one  wavelength 
of  the  undulator  as  one  wavelength  of  light  passes  over  it.  The  dimensionless  coordinate 
z  =  Z/NX  is  normalized  to  the  "slippage  distance"  NX,  the  characteristic  length  along  the 
longitudinal  axis  in  FELs.  It  is  defined  by  the  number  of  optical  wavelengths  that  pass  over  a 
resonant  electron  as  that  electron  traverses  the  undulator  length  L .  Simulations  of  the  FEL  mode 
evolution  take  place  within  a  window  of  width  W  along  z  that  is  an  integral  number  of  slippage 
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distances  long. 

Electrons  in  the  beam  are  governed  by  the  Lorentz  force.  When  the  relativistic  electrons 
(y0> > i )  are  properly  injected  into  the  undulator,  the  transverse  electron  motion  is  periodic  with 
amplitude  and  phase  exp (ik0Z).  The  electron  beam  travels  along  the  undulator  at  an 

average  speed  p0c  where  p0  =  1  -  (1+AT2 )/  2yo .  At  any  dimensionless  time  x  =  p0cr/L ,  electrons  at 
coordinate  z+x  in  the  electron  beam  overlap  the  light  at  site  z  in  the  optical  wave.  The  electron 
phase  relative  to  the  optical  wave  and  undulator  field  is  £  =  (fc+£0)Z-  cor,  and  the  phase  velocity  is 
v  =  £  =  L  [  p2it0  -  *  (1-p2)  ]:  a  resonant  electron  has  v  =  0.  The  Maxwell-Lorentz  equations  [24]  are 

fl(z)  =  -;(z+x)  <  [  1  -  v  (z+z)/2nN  ] 1/2  e-‘5(,+t>  >  (1) 

C  (z+x)  =V  (z+x)  =  8  +  [  1  -  v  (z+x)/2t iN  ]  2  |a(z)  |cos(  C  (z+x)  +  <[>(z ) )  , 

where  the  dimensionless  current  density  is  j(z)  =  8N(neKL)  2p(z)/Yo me2,  p(z)  is  the  actual  particle 
density  at  site  z ,  time  derivatives  are  with  respect  to  the  dimensionless  time  x  (note  that  x  =  0->l 
along  the  undulator  length  L),  and  <...>  is  an  average  over  sample  electrons  at  site  z+x  in  the 
electron  beam.  Since  we  have  assumed  a  slowly-varying  amplitude  and  phase  for  a{z),  the 
sampling  sites  z  can  be  taken  to  be  several  optical  wavelengths  apart.  The  dimensionless  current 
density  j (z)  can  describe  a  short,  peaked  electron  pulse  or  a  long,  flat  beam. 

The  acceleration  term  8  has  been  inserted  to  characterize  a  possible  taper  in  the  undulator 
design  [9].  This  can  be  accomplished  by  tapering  the  wavelength  Xq,  or  the  field  strength  B .  As 
an  example,  a  steady  taper  in  the  undulator  field  strength  gives  8  =  -  AnNK  (1+XT2) LdKtdZ .  No 
taper  in  the  undulator  properties  implies  8  =  0.  The  accelerating  term  8  in  (1)  acts  continuously  to 
restore  the  resonance  condition  as  the  electrons  lose  energy  to  the  laser  light.  The  factors 
[1-v/2rc.V]  in  (1 )  are  close  to  unity  unless  the  electrons  lose  a  significant  fraction  of  their  energy.  A 
large  change  in  v,  comparable  to  2 kN,  can  drastically  change  the  electron-optical  coupling 
strength.  We  take  N  =  30  in  our  examples,  since  recent  oscillator  and  amplifier  experiments  have 
been  in  that  range  [2,3]. 

The  dimensionless  variables  in  (1)  have  been  defined  to  have  meaningful  physical 
interpretations  for  a  wide  range  of  FEL  designs.  Weak  optical  fields  give  values  of  \a  |  <tc,  while 
strong  fields  give  values  of  |  a  |  >n.  The  electron-optical  coupling  is  only  significant  for  |  v  |  not  too 
large,  so  that  the  oscillating  optical  and  undulator  forces  are  nearly  resonant.  The  optical  power 
spectrum  P(k)  is  made  more  relevant  by  expressing  k,  the  optical  wavenumber,  in  terms  of  the 
corresponding  resonance  parameter  v  (k)\  similarly,  the  electron  distribution  function  f(y)  is 


expressed  in  terms  of  v  (7).  The  power  spectrum  F(v  ( k ))  and  the  electron  distribution  function 
/(v(y))  are  more  physically  meaningful  in  terms  of  their  effect  on  the  resonance  condition.  The 
dimensionless  current  density  j  represents  low  single-pass  gain  when  j~  1,  and  high  gain  when 
j»  1.  Typically,  the  FEL  oscillator  configuration  uses  j~1,  while  the  FEL  amplifier  configuration 
uses  y»1.  When  Eqs.  (1)  [24]  are  cast  in  dimensional  form,  they  are  the  same  as  the  Kroll- 
Morton-Rosenbluth  equations  [25]. 

Initial  conditions  are  important  to  the  integration  of  (1).  We  represent  the  beam  at  each  2 
with  a  uniform  distribution  of  initial  phases  £  (x=0)  over  a  2k  range;  the  initial  optical  phase  is 
therefore  arbitrary  and  we  take  <j>(x=0)  =  0  for  each  2.  In  both  the  amplifier  and  oscillator  cases, 
we  start  the  optical  field  at  low  values  a(x=0)  =  a0  ,  well  below  saturation.  The  initial  resonance 
condition  v  (x=0)  =  v0  is  taken  to  be  at  the  maximum  of  the  weak-field  gain  spectrum.  For  small 
current  j  and  weak  fields  a0«  1,  the  gain  found  from  (1)  is  G  =0.135;  at  v0  =  2.6;  for  large 
current  and  weak  fields,  the  gain  is  G  =  exp((J/2)v3^3)/9  at  v0  =  0.  The  FEL  system  would 
normally  develop  significant  coherence  and  power  at  v0  before  going  on  to  saturation. 

3.  Trapped-Particle  Oscillations 

Before  starting  multimode  simulations,  it  is  instructive  to  present  a  phase-space  picture  of 
the  motion  of  trapped  electrons.  Fig.  1  shows  the  ( £,  v )  phase-space  evolution  of  six  sample 
electrons  in  a  strong  optical  field  <z0  =  4k2;  this  field  strength  causes  one  oscillation  of  the  trapped 
electrons.  Many  more  electrons  are  used  in  the  simulation,  but  only  six  are  shown  to  keep  the 
picture  less  cluttered.  At  x  =  0,  the  electron  positions  are  indicated  as  light  grey  spots  started  at 
v0=2.6  and  uniformly  spread  in  phase.  They  are  indicated  by  darker  spots  as  x  increases,  and 
are  black  at  x  =  1.  A  2k  span  in  phase  £  corresponds  to  a  span  along  the  electron  beam  about 
one  optical  wavelength  in  length.  Electron  motion  in  v  corresponds  to  changes  in  electron  energy 
through  v  ~2nN(y-y0)  iy0.  Most  of  the  electrons  shown  in  Fig.  1  are  trapped  in  closed  orbits 
indicating  a  strong  optical  field.  The  gain,  ln(1+G)  where  G  =(|a  \2-a3)lal,  is  plotted  to  the 
right  as  a  function  of  x  along  the  undulator.  As  bunching  develops,  the  gain  and  the  optical  phase 
<j)(x)  grow  as  shown.  The  final  gain  here  is  small,  because  the  FEL  is  well  into  saturation  for  this 
current  density,  j  =5,  and  strong  field  | a  |.  Fig.  1  is  a  single-mode  simulation  since  all  sites  2 
have  been  given  the  same  values  of  v0,  and  0(2)  =  n0. 
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with  the  initial  position  ( C* .  vo  )•  The  case  shown  in  Fig.  1  has  8  =  0  with  small  energy  extraction, 
since  the  changes  in  resonance  are  «2kN  and  N  =  30.  The  synchrotron,  or  trapped-particle, 
oscillation  frequency  is  v5  =  ( a0  sin(  ))1/2  =  a l'2  when  S/a0«  1.  When  the  trapped  electrons 
oscillate  through  a  synchrotron  cycle,  the  gain  in  Fig.  1  also  oscillates  through  one  cycle.  It  is  the 
oscillation  in  the  driving  phase  and  the  FEL  gain  that  causes  the  trapped-particle  frequency  to  be 
imposed  on  the  optical  wave  as  it  slips  over  electrons.  The  sidebands  appear  at  v0  ±  v.s,  so  that 
the  new  FEL  power  is  shifted  from  the  fundamental  wavelength  by  NX!X  =  vs/2kN.  The  shift  has  a 
simple  interpretation;  A XJX  =  "the  number  of  synchrotron  oscillations"/  "the  number  of  undulator 
periods". 

4.  Short  Pulse  Oscillators 

RF  accelerators  inject  a  series  of  short  picosecond  pulses  into  the  FEL  undulator  in  the 
oscillator  configuration  [2],  The  optical  pulses  start  from  spontaneous  emission,  and  bounce 
between  the  mirrors  of  an  optical  resonator  separated  by  a  distance  S  >L .  High-power  saturation 
is  reached  after  several  hundred  passes,  and  the  FEL  works  in  steady-state  for  an  additional  103 
to  1C4  passes.  The  current  density  of  each  short  pulse  j(z)  is  taken  to  be  parabolic  with  the  form 
j{z)~j  ( 1  -2 z2/of)  for  1  z  1  <  cz/^2  and  j{z)~ 0  for  ]z  )>  cr2/V2;  the  FWHM  o2  is  normalized  lo 
the  slippage  distance  NX.  Typically,  RF  accelerators  produce  current  densities  that  give  values  of 
j  in  the  moderate  range  1->  10  and  a2  =  1->  30.  In  addition  to  gain,  there  is  loss  on  each  pass 
due  to  mirror  absorption  and  transmission;  in  the  absence  of  gain,  the  optical  power  decays  as 
e-"IQ  where  n  is  the  pass  number.  Usually,  Q  is  from  10->  200. 

An  important  concept  in  the  short  pulse  FEL  configuration  is  the  matching  of  the  electron 
pulse  repetition  frequency  and  the  bounce  frequency  of  the  light  pulse  in  the  resonator  2 Sc . 
These  frequencies  must  be  closely  synchronized  so  that  each  new  electron  pulse  arrives  at  the 
beginning  of  the  undulato1'  simultaneous  to  the  rebounding  optical  pulse.  Define  d,  the 
"desynchronism,"  as  the  displacement  between  the  pulses  on  each  pass;  if  the  mirrors  are  too 
close  together  by  the  distance  AS,  then  d  =  2A S/NX.  The  desynchronism  d  is  important  to  pulse 
evolution  and  final  saturation.  When  \d\  is  too  large,  the  electron  and  optical  pulses  do  not 
overlap  over  a  sufficient  number  of  passes  and  the  FEL  is  below  threshold  coupling.  Surprisingly, 
when  d  =  0,  exact  synchronism,  the  FEL  is  also  below  threshold  [5-7,11-13,21);  this  effect  has 
been  termed  "laser  lethargy"  [12],  Because  of  slippage,  the  gain  is  preferentially  deposited  on  the 
trailing  edge  of  the  optical  pulse;  therefore,  the  optical  pulse  centroid  actually  travels  slower  than  <. 
in  vacuum,  and  the  frequency  2 Sc  is  overestimated.  To  compensate  for  the  "lethargic"  light,  the 
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path  S  must  be  reduced  by  operating  at  d> 0.  Maximum  power  is  obtained  at  small  d> 0.  The 
final  steady-state  power  P{d)  is  sharply  peaked  at  d  =0,  and  slopes  to  zero  as  d  increases. 


XXX*  FEL  Pulse  Propagation  **** 
J=5  cr2  =10  d=0.02 

Q=10  <5/tT=0  N=30 


Fig.  2.  Several  pictures  show  the  results  of  a  short-pulse  oscillator  simulation  affected  by 
the  trapped-particle  instability.  The  optical  power  P(n)  and  gain  G(n)  evolve  through  «=600 
passes  of  the  oscillator.  The  final  optical  field  amplitude  |a(z)|  and  phase  <) >(z)  are 
modulated  with  the  synchrotron  period  vs  =  2it  while  being  driven  on  each  pass  by  the 
electron  pulse  j(z+ t).  G(v)  is  the  weak-field  gain  spectrum  showing  the  position  of 
maximum  gain  for  the  fundamental;  the  final  optical  power  spectrum  P(v)  has  additional 
sidebands  v0±  vs  due  to  the  instability.  The  final  electron  distribution  /  (v)  is  broad  because 
of  the  strong  fields. 

Figure  2  displays  the  output  of  a  short  pulse  simulation  for  n  =  600  passes  with  j  =  5, 
c,  =  l0,  N  =30,  6  =  0  (no  taper),  c/  =  0.02,  and  Q  =10.  The  Q  is  low  in  this  example  so  that 
sidebands  will  not  mix  further  to  produce  a  chaotic  spectrum.  After  many  passes,  the  initial 


conditions  are  inconsequential,  and  we  look  for  a  final  steady-state  solution  at  saturation.  The 
lower  left  picture  shows  the  electron  pulse  shape  before  (t  =  0)  and  after  (t  =  1)  slippage  of  one 
unit.  The  evolution  of  the  pulses  is  followed  in  a  window  of  width  W  =  16  so  that  -W/2<z<W/2. 
The  upper  left  figures  show  the  final  optical  field  amplitude  \a{z)  |  and  the  optical  phase  profile 
<|>(z).  The  jumps  in  phase  by  2k  are  meaningless,  but  the  slope  of  the  phase  profile  is  a  change  in 
the  laser  light  wavenumber  and  a  resonance  change  Av(£)  =  -<}>'(z)  away  from  v0.  The  grey 
scale  on  §{z)  shows  the  local  frequency  component  matching  the  gain  spectrum  G  (v)  on  the  right. 
The  grey  scale  superimposed  on  the  field  amplitude  above  shows  the  gain  experienced  by  each 
part  of  |a(z)  |  on  the  last  pass  in  the  undulator.  This  grey  scale  matches  the  scale  in  the  gain 
evolution  G  (n)  picture  on  the  right.  The  top  right  figure  follows  the  power  P  (« )  averaged  over  all 
the  sites  in  the  window.  Second  on  the  right  is  G[n)  which  is  the  net  gain  over  the  whole  pulse. 
Third  is  the  final  electron  distribution  /  (v),  and  fifth  is  the  final  optical  power  spectrum  P{v)  found 
from  the  Fourier  transform  of  a{z).  The  fourth  picture  G(v)  is  the  weak-field,  single-mode  gain 
spectpjm  for  the  current  density  j  on  a  single  pass;  it  is  included  for  reference. 

The  output  of  the  simulation  shows  that  the  FEL  oscillator  has  reached  saturation  with  peak 
fields  near  \a  |  =  40.  On  each  pass,  those  electrons  that  are  near  the  center  of  the  pulse,  and 
have  initial  phases  that  allow  them  to  become  trapped,  encounter  fields  strong  enough  to  cause 
about  one  synchrotron  oscillation.  Normal  saturation  takes  place  early  and  the  power  is  steady  for 
about  350  passes;  then  the  power  increases  again  with  the  onset  of  the  instability.  The  extra 
power  is  in  the  sidebands  and  the  fundamental  remains  saturated  [17,23].  Strong  optical  fields 
and  the  resulting  trapped-particle  oscillations  must  occur  before  there  is  significant  gain  at  the 
sideband  frequency  vs.  A  simple  calculation  shows  that  the  field  component  at  the  sideband 
frequency  vs  grows  as 

a  (vs)  =  as  exp  (0'5/4v5)  1/2  t)  (3) 

where  vs  =  \a  j 1/2  at  saturation,  as  is  the  initial  sideband  field  at  vs,  and  js  is  the  trapped  fraction 
of  j\  typically,  js  =  (0.1 ->0.5);'.  In  this  simulation,  OV4vs)v2  =  0-3  and  the  loss  rate  is 
-  (2  Q  T1  =-0.05,  so  the  sideband  gain  is  above  threshold.  The  sideband  structure  is  clear  in 
|a(z)  | ,  and  the  power  spectrum  P(v).  The  modulation  length  is  close  to  the  slippage  distance 
indicating  one  synchrotron  oscillation  of  the  trapped  electrons  each  pass,  and  the  sideband  P(v) 
occurs  at  vs  =  2k  above  the  fundamental.  At  the  ends  of  the  pulse  where  the  fields  are  weaker, 
the  sideband  period  is  slightly  longer,  in  agreement  with  (2).  The  full-width  of  the  electron 
distribution  is  given  by  the  height  of  the  closed-orbit  region  in  phase-space,  4 1  a  \ 1/2 


The  general  features  of  the  trapped-particle  instability  in  short-pulse  FELs  are  briefly 

outlined  below  [21];  many  features  have  now  been  observed  in  experiments  [2,8,1 1], 

1.  At  small  d> 0,  the  FEL  usually  reaches  power  levels  large  enough  to  cause  the  trapped- 
particle  instability.  In  this  case,  the  optical  pulse  will  be  centered  on  the  electron  pulse,  and 
will  have  sharp  spikes  due  to  the  instability.  This  gives  a  broad,  possibly  chaotic,  optical 
power  spectrum,  and  a  broad  electron  distribution. 

2.  At  large  d,  the  steady-state  power  is  smaller  due  to  the  reduced  coupling,  and  the  trapped- 
particle  instability  is  less  likely  to  occur.  The  final  optical  power  spectrum  is  narrow  and  in  a 
single-mode;  the  final  electron  distribution  is  narrow  due  to  the  weak  optical  fields.  Since  the 
optical  pulse  is  advanced  by  a  large  d  on  every  pass,  the  center  of  the  optical  pulse  may 
actually  be  ahead  of  the  electron  pulse.  If  the  trapped-particle  instability  is  undesirable  in  an 
FEL  application,  it  may  be  easily  removed  by  increasing  d. 

3.  When  d  is  in  the  intermediate  range,  we  have  often  observed  limit-cycle  behavior  in  the 
simulations  [5,6,13,21],  In  this  case,  the  pulse  continually  changes  shape  while  the 
trapped-particle  instability  creates  new  subpulses. 

4.  Increasing  the  current  density  j  or  the  resonator  Q  increases  the  steady-state  power;  this 
increases  the  synchrotron  frequency  and  the  sideband  gain.  The  addition  of  sideband 
power  is  cumulative,  since  the  presence  of  a  strong  sideband  again  increases  the  steady- 
state  power.  Simple  modulation  as  shown  in  Fig.  2  is  not  the  usual  case;  most  parameter 
choices  lead  to  a  chaotic  optical  spectrum  with  many  random  spikes  due  to  further  sideband 
mixing  [22], 

5.  If  taper  is  introduced  into  the  undulator  design,  the  synchrotron  frequency  is  only  slightly 
modified.  But,  the  depth  of  the  strong-field  potential  wells  is  significantly  reduced,  and  a 
smaller  fraction  of  the  electrons  are  trapped  in  closed  harmonic  orbits.  Tapered  undulators 
reduce  the  sideband  gain  and  tend  to  suppress  the  trapped-particle  instability  [2,13,22,23]. 

6.  When  the  pulse  length  a,  is  near  unity,  the  short-pulse  effects  can  be  dramatic  [11]. 
Adjusting  d  affects  the  total  power,  the  electron  distribution,  and  the  optical  pulse  structure 
as  outlined  above.  Since  the  sideband  modulation  first  appears  near  the  slippage  distance 
(normalized  to  unity),  a  short  pulse  may  not  be  long  enough  to  support  modulation  at  the 
synchrotron  frequency.  Depending  on  d  the  optical  pulse  can  be  made  significantly  shorter, 
or  longer,  than  the  picosecond  electron  pulse. 
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5.  Long-Pulse  Oscillators,  "Wrapped-YVindows" 

In  long  electron  pulses  (  ct,»1  )  distant  slippage  sections  of  length  Az=1  are  largely 
uncoupled  and  evolve  independently  in  a  window  tV>oJ»1.  Simulations  can  become  time 
consuming  and  require  large  amounts  of  memory.  The  straight-forward  method  used  in  Fig.  2 
becomes  wasteful  because  many  of  the  slippage  sections  have  similar  evolutions.  A  better  way  to 
solve  the  problem  [17,22,23]  is  by  sampling  a  smaller  window  ^<<0,  with  periodic  bounaary 
conditions  ( (z-VV/2)  =  C,  (z+WI2).  This  restricts  the  number  of  modes  that  can  be  examined. 
Using  a  number  of  sites  Nw  in  the  window  W  we  only  represent  modes  v,  =  v0  -  (2k/  W)(  1-NWI2 ) 

where  the  integer  /=  0,1,2 . Wiv-1;  lhe  mode  spacing  is  Av,  =  2 n/W.  With  the  periodic 

boundary  conditions  applied  to  the  ends  of  the  windows,  we  call  them  "wrapped-window" 
simulations. 

Fig.  3  shows  the  result  of  a  simulation  in  a  window  IV  =  4.  The  FEL  is  described  by  j  =  5, 
Q  =  10,  N  =  30,  and  8  =  0  as  in  Fig.  2.  Without  pulses,  d  and  cr,  don’t  enter  the  problem.  The 
individual  pictures  are  the  same  as  in  Fig.  2  except  for  the  reference  to  the  pulse  shape  and 
slippage.  The  electron  current  density  j(z)  is  constant  along  the  window  with  a  small  amount  of 
random  noise  Sj  =  0.01  j  at  each  site.  Without  noise,  every  site  z  would  evolve  identically  and  no 
spectral  features  would  develop,  but  specific  sources  of  noise  can  vary  from  one  experiment  to  the 
next.  The  dimensionless  current  density  j  depends  on  several  physical  quantities  which  could 
contribute  noise.  Fluctuations  of  the  electron  particle  density  p(z),  or  undulator  errors  in  K2,  are 
both  candidates. 

After  n  =  600  passes,  the  steady-state  optical  fields  peak  near  |  a  |  =  40  just  as  in  Fig.  2,  and 
the  final  power  spectrum  P{\)  is  nearly  identical  to  Fig.  2.  Yet,  this  wrapped-window  simulation 
takes  only  1/4  the  computer  time  of  the  pulse  simulation;  a  smaller  window  IV  =  2  would  save 
more  time  and  give  the  same  result. 
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Fig.  3.  With  periodic  boundary  conditions  in  a  "wrapped-window"  simulation,  the  modulation 
and  spectral  lines  resulting  from  the  trapped-particle  instability  are  made  more  clear.  The 
complicating  features  of  short  pulses,  seen  in  Fig.  2,  are  not  present. 
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The  general  features  of  wrapped-window  simulations  are  the  same  as  the  short-pulse  cases 
except  for  desynchronism.  An  advantage  of  the  wrapped-windows  is  the  reduced  number  of 
dimensionless  parameters. 

1.  As  with  pulse  simulations,  increases  in  j  or  Q  make  the  trapped-particle  instability  more 
severe.  It  is  unusual  to  find  such  even  modulation  as  shown  in  Fig.  3.  The  instability 
increases  output  power  by  increasing  the  laser  bandwidth. 
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2.  When  there  is  taper  in  the  undulator  design,  the  instability  is  less  severe  for  the  same 
reasons  stated  above. 

3.  Some  source  of  noise  is  essential,  but  details  are  not  important  in  the  oscillator  case.  For 
pulse  simulations  no  random  noise  is  introduced,  since  the  spectral  features  of  the  electron 
pulse  structure  are  sufficient.  In  fact,  we  typically  see  the  sideband  modulations  start  at  the 
ends  of  the  electron  pulse  where  the  Fourier  components  are  larger.  In  the  FEL  oscillator 
case,  there  are  one  or  more  synchrotron  oscillations  each  pass.  Over  many  hundreds  of 
passes,  the  stored  optical  wave  "sees"  many  synchrotron  oscillations,  and  even  modest 
sideband  gain  above  threshold  gives  enormous  growth  from  a  small  amount  of  noise.  The 
resulting  steady-state  features  are  not  affected  by  the  details  of  the  noise  source  employed. 

6.  High-Current  FEL  Amplifiers 

The  FEL  amplifier  uses  a  high  current  density  j»  1  to  reach  high-power  in  a  single  pass.  It 
is  possible  to  attain  values  of  j  =  1 04— >1 05  with  the  use  of  induction  linac  accelerators  and  long 
undulators  [3].  The  electron  pulses  are  much  longer  than  the  slippage  distance,  a,  ->  °°,  and  the 
wrapped-window  simulation  method  becomes  essential.  The  FEL  growth  rates  are  so  large  that 
electrons  are  trapped  early  in  the  undulator  and  begin  executing  synchrotron  oscillations.  There 
can  be  from  several  to  many  tens  of  synchrotron  oscillations  along  the  undulator,  but  these  are  the 
only  synchrotron  oscillations  experienced  by  the  optical  field.  Even  with  far  fewer  synchrotron 
cycles  than  in  the  oscillator  case,  the  large  current  density  j  can  give  significant  sideband  gain 
once  trapping  occurs. 

Fig.  4  shows  the  final  simulation  output  for  an  FEL  with  high  current  j  =5x104,  no  taper 
8=  0,  an  initial  field  a0  =  20  starting  at  v0  =  0,  and  N  =  30.  A  small  amount  of  noise  is  introduced 
over  a  wide  range  of  frequencies  using  8y  =  0.01  j ;  in  addition,  there  is  a  specific  and  comparable 
contribution  from  an  initial  sideband  field  as  =  0.1  at  vs  =  16tc.  A  high-current  beam  with  some  of 
the  electrons  trapped  by  strong  fields  will  provide  a  source  for  as;  small  fluctuations  from  shot 
noise  or  spontaneous  emission  in  the  average  <...>  lead  to  significant  contributions  to  as  because 
j  is  large.  The  power  evolution  /’(x)  shows  about  eight  synchrotron  oscillations;  these  oscillations 
in  a(x)  are  converted  to  the  a[z)  modulation  through  the  trapped-particle  instability.  The  sideband 
appears  to  grow  at  a  rate  consistent  with  (3)  as  soon  as  the  field  |  a  (x)  |  becomes  sufficiently 
strong  to  cause  trapped  particle  oscillations.  The  final  field  modulation  and  sideband  in  the  power 
spectrum  are  clear. 
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Fig.  4.  In  the  high-current  FEL  amplifier,  the  trapped-particle  instability  must  occur  in  one 
pass.  In  much  stronger  optical  fields,  the  trapped  electrons  experience  eight  synchrotron 
oscillations,  and  impose  a  corresponding  modulation  onto  the  optical  wave  envelope. 
Sidebands  are  spaced  at  vs  =  1  6k. 

The  general  characteristics  of  wrapped-window  simulations  for  FEL  amplifiers  are  different 
in  many  ways  from  the  oscillator  case.  There  is  no  resonator  Q  to  consider,  and  sources  of  noise 
are  much  more  important  to  FEL  amplifiers  than  to  FEL  oscillators. 

1.  As  with  the  oscillator  cases,  large  j  makes  the  trapped-particle  instability  more  severe. 
However,  it  is  only  the  large  values  of  j  that  have  a  possibility  of  developing  a  sideband 
comparable  in  power  to  the  fundamental.  Using  (3)  we  can  evaluate  the  current  j  needed 
for  the  sideband  field  as  to  become  equal  to  the  fundamental  amplitude  a0  =  v|.  We 
estimate  that  most  of  the  current  is  trapped  in  harmonic  orbits,  js  ~  0.75 j  ,  in  the  high 


current  case.  With  high  current,  the  potential  well  that  traps  electrons  grows  rapidly  in  height 
and  forces  most  of  the  current  into  harmonic  orbits. 

j*  ~  20  Vy  ln2(  vs  A/«7)  (4) 

Assume  that  at  least  a  few  synchrotron  oscillations  (  say  v5  =  1 67c )  are  necessary  for  the 
optical  wave  to  pick  up  the  synchrotron  frequency.  The  initial  sideband  field  is  estimated  at 
=  0.1  for  these  high  current  cases;  j*  only  has  a  logarithmic  dependence  on  as  so  that 
the  estimate  is  not  crucial.  Then,  we  find  that  a  characteristic  current  density  j*  =  3x1 04  is 
needed  before  sidebands  can  significantly  alter  the  final  FEL  amplifier  spectrum,  and 
currents  below  j *  are  not  as  susceptible  to  the  trapped-particle  instability.  The  example  in 
Fig.  4  uses  j  =  5x104  >  j*  for  this  reason.  We  found  the  trapped-particle  instability  difficult 
to  excite  when 

2.  When  there  is  taper  in  the  undulator  design,  the  instability  is  less  severe  for  the  same 
reasons  as  with  the  FEL  oscillators.  Less  current  is  trapped  and  the  synchrotron  oscillations 
have  a  smaller  amplitude. 

3.  Unlike  the  FEL  oscillators,  the  spectral  features  of  noise  are  important  to  sideband  growth  in 
the  FEL  amplifier.  Suppression  of  input  noise  at  the  sideband  frequency  could  be  an 
important  method  for  avoiding  the  trapped-particle  instability  in  amplifiers.  However,  shot 
noise  and  spontaneous  emission  cannot  be  eliminated,  and  are  significant  contributions  at 
these  current  densities. 

4.  The  high  current  necessary  for  the  trapped-particle  instability,  also  plays  a  role  in 
suppressing  it.  For  such  large  j  >j * ,  there  is  no  steady-state  achieved  even  in  strong  fields 
|  a  | ;  the  power  continues  to  increase  and  so  does  the  synchrotron  frequency.  This  can  be 
seen  in  P(x)  of  Fig.  4.  As  the  synchrotron  frequency  changes,  the  sideband  power  at  a 
particular  frequency  does  not  accumulate. 

In  conclusion,  we  note  that  the  trapped-particle  instability  is  more  likely  to  be  found  in  the 
FEL  oscillator  than  in  the  FEL  amplifier,  even  though  it  was  first  predicted  for  the  amplifier  case. 
This  is  due  to  the  large  number  of  synchrotron  oscillations  experienced  by  the  light  stored  in  the 
oscillator  resonator  for  many  passes.  However,  the  adjustable  desynchronism  d  (and  even  a 
wavelength  selective  g{v(X))  can  make  it  relatively  easy  to  reduce  the  trapped-particle  instability 
in  the  oscillator  case,  while  the  FEL  amplifier  allows  no  such  control  in  just  one  pass.  To  suppress 
sidebands  in  the  amplifier,  it  may  prove  most  fruitful  to  explore  particular  taper  designs  5(Z),  the 


input  wavelength  v0(X)  and  field  strength  a0,  and  controllable  sources  of  noise  like  undulator 
errors  and  current  fluctuations.  It  should  also  be  remembered  that  while  the  trapped-particle 
instability  has  already  been  observed  in  the  FEL  oscillator,  it  has  not  yet  been  found  in  the  FEL 
amplifier. 
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ABSTRACT 


In  high-power  free-electron  lasers,  the  electrons  can  become  trapped  in  deep 
potential  wells  formed  by  the  combined  opv,cal  and  undulator  field  forces.  The 
trapped  current  oscillates  at  the  synchrotron  frequency,  and  can  drive  the  optical 
wave  at  sideband  frequencies  around  the  fundamental.  This  "trapped-particle 
instability"  can  occur  in  both  the  oscillator  and  amplifier  configurations. 


1.  Introduction 

Many  free-electron  laser  (FEL)  oscillators  [1],  and  FEL  amplifiers  [2]  are  designed  to  pro¬ 
duce  large  optical  power.  In  both  configurations,  a  relativistic,  high-current  electron  beam 
amplifies  a  copropagating  electromagnetic  wave  as  they  pass  through  a  periodic,  transverse  mag¬ 
netic  field  undulator.  The  trapped-particle  instability  starts  when  electrons  become  trapped  in 
deep  potential  wells  formed  by  the  combined  optical  and  undulator  fields.  Electrons  near  the  bot¬ 
tom  of  the  well  oscillate  in  harmonic  orbits  at  the  synchrotron  frequency  causing  the  carrier  wave 
to  develop  sidebands.  The  "trapped-particle  instability"  was  first  predicted  using  a  single-mode 
theory  describing  high-power,  low-gain  FEL  oscillators  with  a  tapered  undulator  [3];  the  tapered 
undulator  design  reaches  high  power  by  trapping  electrons  in  potential  wells  that  remain  resonant 
as  the  electrons  lose  energy  [4],  At  the  same  time,  multimode  simulations  of  shori-pulse  FEL 
oscillators  observed  the  effects  of  the  instability  as  pulse  modulation  [5-7],  Recently,  the  instability 
has  been  observed  in  the  high-power  FEL  oscillator  at  LANL  [1],  and  possibly  in  the 
TRW/Stanford  oscillator  [8],  It  has  been  termed  the  sideband,  Kroll-Rosenbluth,  synchrotron, 
Raman,  and  the  trapped-particle  instability  [9-23].  Here,  the  multimode  sideband  theory  is 
reviewed  with  examples  from  short-pulse  FELs,  oscillators,  and  amplifiers.  Common  features  and 
differences  are  discussed. 
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2.  Multimode  Simulation  Theory 


The  electrons  injected  into  the  FEL  undulator  evolve  due  to  the  combined  optical  and  undu- 
lator  electromagnetic  fields.  When  the  relativistic  electrons  of  energy  y< <?ncz  (y0>>1)  are  properly 
injected  into  the  undulator,  the  transverse  electron  motion  is  periodic  with  amplitude  K X$  /2rcy0  and 
phase  exp(iJk0Z)  where  K  =  eBX0l2nmczl  e  is  the  electron  charge  magnitude,  B  is  the  rms  undu¬ 
lator  field  strength,  Xc  =  2n/k0  is  the  undulator  wavelength,  m  is  the  etectron  mass,  and  c  is  the 
speed  of  light  in  vacuum.  The  average  speed  of  the  beam  along  the  Z-axis  is  p0c  where 
p0  =  1  -  (1+A-2 )/  2y<j.  We  follow  the  beam  evolution  with  the  dimensionless  time  x  =  p 0ct/L  =  0->1 
along  the  undulator  length L  =  NXc,m\hN  periods. 

The  light  wave  evolves  in  the  presence  of  the  electron  beam  according  to  the  transverse 
wave  equation.  The  optical  carrier  wave  with  frequency  co  =  J tc  =2 ncfX  has  the  single-mode 
phase  exp[i(/tZ-cof)]  and  a  complex  slowly-varying  coefficient  a(z)  =  aK(z)  +  ia,[z)-  |a(z)|e,i:,(,) 
evaluated  at  many  discrete  sites  z.  The  dimensionless  coordinate  z  is  the  ratio  Z/NX  where  the 
"slippage  distance"  NX  is  defined  by  the  number  of  optical  wavelengths  that  pass  over  a  resonant 
electron  as  that  electron  traverses  the  undulator  length  L .  With  the  slowly-varying  amplitude  and 
phase  approximation,  and  the  coordinate  change  Z-xZ+ct ,  the  wave  operator  reduces  to  a  single 
derivative  in  time.  The  dimensionless  optical  field  envelope  is  a{z)  =  4nNeKLE{z)/y$mcz  where 
E{z)  is  the  complex  optical  electric  field,  and  y0mc2  is  the  resonant  electron  energy, 
yo  =  k(1+Kz)/2k0.  Simulations  take  place  within  a  window  of  width  W  along  z  that  is  an  integral 
number  of  slippage  distances  long.  At  any  x,  electrons  at  coordinate  z  +  x  in  the  electron  beam 
overlap  the  light  at  coordinate  z  in  the  optical  wave.  The  electron  phase  relative  to  the  optical 
wave  and  undulator  fields  is  l^  =  [k+k0)Z-  cor ,  and  the  phase  velocity  is  v  =  £  =  L[  (/.-+£  0)P,  -k  ]. 


The  self-consistent  electron  and  optical  equations  [24]  are 

v(z+x)  =  [  1  -  v  (z+x)/2nW  ] 2  ( a*(z)cos(  C  (z+x) )  -  a/(z)sin(  £  (z+x) ) )  (1) 

£(z+t)  =  v(z+x)  (2) 

ri*(z)  =- j'(z+x)  <  [  1  -v  (z+x)/2irN  ]  1/2cos(  £  (z+x) )  >  (3) 

«/(z)  =y  (z+x)  <  [  1  -  v  (z+x)/2tc/V  ]  v2sin(  C  (z+"t) )  >  (4) 


where  the  dimensionless  current  density  is  j[z)  =  8N (c KKL)zp{z)lyzmc2 ,  p(z)  is  the  actual  particle 
density  at  site  z,  time  derivatives  are  with  respect  to  the  dimensionless  time  x,  and  <...>  is  an 
average  over  sample  electrons  at  site  z  +  x  in  the  electron  beam. 


In  order  to  maintain  strong  coupling  between  the  electrons  and  light,  the  phase  velocity  v 
cannot  be  too  far  from  resonance,  v  =  0.  The  initial  resonance  condition  v  (t=0)  =  v0  is  taken  to  be 
at  the  maximum  of  the  weak-field  gain  spectrum.  In  both  the  amplifier  and  oscillator  cases,  we 
start  the  optical  field  at  low  values  a(x=0)  =  ao,  well  below  saturation.  The  dimensionless  current 
density  j  represents  low  single-pass  gain  when  j~  1,  and  high  gain  when  ;'» 1 .  Typically,  the  FEU 
oscillator  configuration  uses  j~  1,  while  the  FEL  amplifier  configuration  uses  y»1.  Weak  optical 
fields  give  values  of  \a  |  <rc,  while  strong  fields  give  values  of  \a  \  >k.  For  small  current  and  weak 
fields,  maximum  gain  is  at  v0  =  2.6;  for  large  current  and  weak  fields,  maximum  gain  is  at  v0  =  0. 
The  electron  beam  at  each  z  has  a  uniform  distribution  of  initial  phases  £  (t=0)  over  a  2k  range  so 
that  the  initial  optical  phase  is  arbitrary;  we  take  4>(x=0)  =  0  at  each  z.  The  factors  [1-v/27iW]  in 
(l)-(4;  are  close  to  unity  unless  the  electrons  lose  a  significant  fraction  of  their  energy  and 
become  less  relativistic.  Electron  motion  in  terms  of  v  corresponds  to  changes  in  the  electron 
energy  through  the  relation  v  =  4kN ( y- y0 )  /  y0.  We  take  N  =  50  in  our  examples  as  a  common 
value  typical  of  both  amplifiers  and  oscillators.  The  final  optical  power  spectrum  P(k)  is  made 
more  relevant  by  expressing  k ,  the  optical  wavenumber,  in  terms  of  the  corresponding  resonance 
parameter  v  [k)\  similarly,  the  electron  distribution  function /  (y)  is  expressed  in  terms  of  v  (y).  The 
power  spectrum  P(v  ( k ))  and  the  electron  distribution  function  /  (v  (y))  are  more  physically  mean¬ 
ingful  in  terms  of  their  affect  on  the  resonance  condition. 

The  single-mode  version  of  (1)-(4)  is  obtained  by  removing  all  z  dependence;  all  sites  are 
given  the  same  initial  values  of  £0.  v0,  and  a{z)  =  a0.  An  estimate  of  the  trapped-particle  motion  is 
obtained  by  considering  the  electrons  in  harmonic  orbits  near  the  stable  fixed-point  =  n/2.  With 
j  not  too  large,  and  and  small  energy  extraction  (N»  1),  | a  |  remains  approximately  constant 
during  the  synchrotron  oscillations  at  saturation.  Then,  the  motion  of  a  trapped  electron  is 
CM  =  C+  ( vo  ,vs)  sin(  vst  )  with  the  initial  position  (  C* .  v0)-  The  synchrotron  or  trapped-particle 
oscillation  frequency  is  \s  =  ao2.  When  the  trapped  electrons  oscillate  through  a  synchrotron 
cycle,  part  of  the  current  driving  the  optical  field  in  (3)  and  (4)  also  oscillates  through  one  cycle.  It 
is  the  oscillation  of  the  driving  phase  in  the  average  <...>  that  causes  the  trapped-particle  fre¬ 
quency  to  be  imposed  on  the  optical  wave  as  it  slips  over  electrons.  The  sidebands  appear  at 
v0  ±  vs ,  so  that  the  new  FEL  power  is  shifted  from  the  fundamental  wavelength  by  AXfk  =  vs/2kN  . 
The  shift  has  a  simple  interpretation;  /\X/X=  "the  number  of  synchrotron  oscillations"/  "the  number 
of  undulator  periods". 


3.  The  Trapped-Particle  Instability  in  Short-Pulse  Oscillators 

An  FEL  oscillator  that  is  powered  by  an  RF  accelerator  injects  a  series  of  short  picosecond 
electron  pulses  into  the  undulator  while  the  optical  pulse  bounces  between  mirrors  separated  by  a 
distance  S>L.  High-power  saturation  is  reached  after  many  passes,  and  the  FEL  continues  to 
work  for  an  additional  103  to  104  passes.  The  current  density  of  each  short  pulse  j(z)  is  taken  to 
be  parabolic  with  the  form  j{z)  =j  ( 1  - 2 z2/ct22 )  for  |z  |  <  a2/VZ  and  j(z)= 0  for  |z  | >  az/V?;  the 
length  ct2  is  normalized  by  the  slippage  distance  NX.  Typically,  RF  accelerators  produce  current 
densities  which  give  values  of  j  in  the  moderate  range  1  — >  100  and  az  =  1->  30.  The  loss  on  each 
pass  (due  to  mirror  absorption  and  transmission)  is  described  by  e~n/Q  where  n  is  the  pass 
number.  Usually,  Q  is  from  2->  200. 

The  repetition  frequency  of  successive  electron  pulses  must  be  matched  to  the  bounce  fre¬ 
quency  of  the  light  pulse,  2Sc.  When  synchronized,  each  new  electron  pulse  arrives  at  the  begin¬ 
ning  of  the  undulator  simultaneous  to  the  rebounding  optical  pulse.  The  "desynchronism," 
d  =  2AS'NX,  is  the  displacement  between  the  pulses  after  each  pass  when  the  mirrors  are 
separated  by  S  -  AS.  If  \d\  is  too  large,  the  electron  and  optical  pulses  do  not  overlap  for  a 
sufficient  number  of  passes  and  the  FEL  operates  below  threshold  coupling.  If  d  -  0,  exact  syn¬ 
chronism,  the  FEL  is  also  below  threshold  [5-7,11-14,21]  due  to  an  effect  termed  "laser  lethargy" 
[12].  Because  of  slippage,  gain  is  preferentially  deposited  on  the  trailing  edge  of  the  optical  pulse 
causing  the  optical  pulse  centroid  to  travel  slower  than  c  in  vacuum;  therefore,  the  optical  bounce 
frequency  25c  is  overestimated.  To  compensate  for  the  "lethargic"  light,  the  path  5  must  be 
reduced  by  operating  at  d> 0. 


FEL  Pulse  Propagation 
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Fig.  1.  Driven  by  the  short  electron  pulse  j[z+i),  the  trapped-particle  instability  has  modulated 
the  optical  pulse  envelope  a{z)  at  the  synchrotron  period  vs  =  2k.  The  total  optical  power 
P{n)  and  gain  G(n)  reach  steady-state  after  «=400  passes  through  the  oscillator.  G(v)  is 
the  weak-field  gain  spectrum,  included  for  reference.  The  final  optical  power  spectrum  P(v) 
has  an  additional  sideband  at  v0  +  v,s  due  to  the  instability.  The  final  electron  distribution 
/(v)  is  broad  because  of  the  strong  fields. 

A  short-pulse  simulation  is  shown  in  Fig.  1.  Steady-state  is  reached  after  «  =400  passes 
with  moderately  high  current  j  =  20,  an  intermediate  pulse  length  a,  =  12,  N  =  50,  desynchronism 
d  =  0.04,  and  0=3.  Large  output  coupling  is  used  ( low  Q  )  because  the  gain  is  far  above  thres¬ 
hold.  The  top  right  figure  follows  the  power  P(n)  averaged  over  the  whole  window.  Second  on 
the  right  is  the  net  gain  G{n)  experienced  by  the  pulse  each  pass.  The  lower  left  picture  shows 
the  electron  pulse  shape  before  (x  =  0)  and  after  (x  =  1)  the  slippage  of  one  unit  per  pass  in  a  win¬ 
dow  of  width  W  =  20.  The  upper  left  figures  show  the  final  optical  field  amplitude  |  a  (z )  |  and  the 
optical  phase  profile  <j>(z)  plotted  along  z  where  -W /2<z  <W  12.  The  jumps  in  phase  of  2n  are 
meaningless,  but  the  slope  of  the  phase  profile  indicates  a  local  change  in  the  laser  light 
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wavenumber  and  a  resonance  change  Av(ifc)  =  -(S>'(2)  away  from  v0.  The  grey  scale  on  <j>(z) 
shows  the  local  frequency  component  matching  the  gain  spectrum  G  (v)  on  the  right.  The  grey 
scale  superimposed  on  the  field  amplitude  shows  the  local  gain  experienced  by  each  part  of  a(z) 
on  the  last  pass  in  the  undulator.  This  grey  scale  matches  the  scale  in  G  (n)  pictured  on  the  right. 
Third  on  the  right  is  the  final  electron  distribution  /  (v),  and  fifth  is  the  final  optical  power  spectrum 
P  (v)  found  from  the  Fourier  transform  of  a  (z).  The  fourth  picture  G  (v)  is  the  weak-field  gain  spec¬ 
trum  for  the  current  density  j  =20. 

The  simulation  has  reached  saturation  with  peak  fields  near  \a  |  =60.  Each  pass,  those 
electrons  near  the  center  of  the  pulse,  that  become  trapped,  encounter  fields  strong  enough  to 
cause  about  one  synchrotron  oscillation  since  -fa  =  2n.  Normal  saturation  takes  place  early  and 
the  power  is  steady  for  about  50  passes;  then  the  power  increases  again  as  the  sideband  grows. 
Strong  optical  fields  and  the  resulting  trapped-particle  oscillations  must  occur  before  there  is 
significant  gain  at  the  sideband  frequency  vs.  The  extra  power  is  in  the  sidebands  and  the  funda¬ 
mental  remains  saturated  [17,23],  A  simple  calculation  shows  that  the  field  component  at  the 
sideband  frequency  grows  as  a  (t)  =  as  exp  ((y/4vs) 1/2 z)  where  vs  =  \a  \ 1/2  at  saturation,  and  as  is 
the  initial  sideband  field  at  vs.  In  this  simulation,  (;74vs)1/2  =  0.8  and  the  loss  rate  is 
-(2  2)-1  =-0.2,  so  the  sideband  gain  is  above  threshold.  The  sideband  structure  is  clear  in 
|a(z)|,  and  the  power  spectrum  P(v).  The  modulation  length  is  close  to  the  slippage  distance 
indicating  v5  =2jc,  and  the  sideband  in  P(v)  occurs  at  Av  =  2ti  above  the  fundamental.  The  full- 
width  of  the  electron  distribution  is  given  by  the  height  of  the  closed-orbit  region  in  phase-space, 
4 |a  I1'2. 

The  characteristics  of  the  trapped-particle  instability  in  short-pulse  FELs  are  briefly  outlined 
below  [21].  Most  characteristics  have  now  been  confirmed  by  experiments  [1 ,8,1 1]. 

1.  At  small  d> 0,  the  FEL  usually  reaches  power  levels  large  enough  to  cause  the  trapped- 
particle  instability.  This  gives  a  broad,  possibly  chaotic,  optical  power  spectrum,  and  a 
broad  electron  distribution. 

2.  At  large  d,  the  steady-state  power  is  small  due  to  the  reduced  coupling,  and  the  trapped- 
particle  instability  is  less  likely  to  occur.  The  final  optical  power  spectrum  is  narrow  and  in  a 
single-mode;  the  final  electron  distribution  is  narrow  due  to  the  weak  optical  fields. 


3.  When  d  is  in  the  intermediate  range,  we  have  often  observed  limit-cycle  behavior  in  the 
simulations  [5,6,13,21].  In  this  case,  the  pulse  continually  changes  shape  while  the 
trapped-particle  instability  creates  new  subpulses. 

4.  Increasing  the  current  density  j  or  the  resonator  Q  increases  the  steady-state  power,  the 
synchrotron  frequency,  and  the  sideband  gain.  The  addition  of  sideband  power  is  cumula¬ 
tive,  since  the  presence  of  a  strong  sideband  again  increases  the  steady-state  power. 

5.  If  taper  is  introduce  into  the  undulator  design,  the  synchrotron  frequency  is  only  slightly 
modified,  and  the  sideband  gain  is  reduced  [13,22,23]. 

6.  When  the  pulse  length  a,  is  near  unity,  the  short-pulse  effects  can  be  dramatic  [11].  Since 
the  sideband  modulation  first  appears  near  the  slippage  distance,  a  short  pulse  may  not  be 
long  enough  to  support  modulation  at  the  synchrotron  frequency.  The  optical  pulse  can  be 
made  significantly  shorter  or  longer  than  the  picosecond  electron  pulse  by  adjusting  d . 

4.  Simulations  in  a  "  Wrapped-Window" 

If  the  electron  pulse  is  not  short  (a,»1)  it  is  prudent  to  simulate  the  FEL  by  sampling  a 
smaller  window  W«a,  with  periodic  boundary  conditions  such  that  C  (z-W/2)  =  £  (z+W/2);  these 
can  be  called  "wrapped-window"  simulations  [17,22,23].  Using  a  number  of  sites  Nw  in  a  window 
of  width  W  we  follow  a  restricted  number  of  modes  v,  =v0-  [2kIW){  1-NWI2)  where 
/  =  0, 1 , 2, ....  Nw- 1 ;  the  mode  spacing  is  Av  =  2n/W.  The  desynchronism  d  and  the  pulse  length 
don’t  enter  the  problem,  and  the  current  density  j{z)=j  is  constant  along  the  window. 

The  FEL  simulation  in  Fig.  2  uses  j  =  20,  2=3,  and  N  =50  (as  in  Fig.  1),  but  for  a  long 
pulse  ct,»W'=  4.  The  individual  pictures  are  the  same  as  in  Fig.  1  except  for  the  reference  to 
pulse  shape.  Without  some  noise  source,  no  power  would  develop  at  frequencies  other  than  the 
fundamental  and  every  site  z  would  evolve  identically.  Specific  sources  of  noise  can  vary  from 
one  experiment  to  the  next.  In  Fig.  2,  the  initial  electron  phases  are  uniformly  spread  over  a  27t 
range,  but  with  an  additional  random  phase  of  zero  mean  and  standard  deviation  5C  =  I.OxlCT6. 
Electron  shot  noise  is  a  typical  source  of  this  kind  of  noise,  In  the  pulse  simulations  no  random 
noise  is  introduced,  since  the  spectral  features  of  ;'(z)  are  sufficient.  After  n  =400  passes,  the 
steady-state  optical  fields  peak  at  \a  |  =  60  just  as  in  Fig.  1 ,  and  the  final  power  spectrum  P  (v)  is 
similar  to  Fig.  1,  but  without  the  short  pulse  features. 
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Fig.  2.  The  "wrapped-window"  simulation  gives  a  clear  example  of  the  trapped-particle  instability 
in  the  FEL  oscillator  without  the  complicating  features  of  short  pulses. 

Wrapped-window  simulations  have  many  of  the  same  characteristics  as  for  pulses,  but 
without  the  effects  of  desynchronism.  As  cr,-»  «>  the  wrapped-window  approach  is  the  only  viable 
solution  to  the  FEL  multimode  simulation. 

1.  Increases  in  j  or  Q  make  the  trapped-particle  instability  more  severe.  It  is  unusual  to  find 
the  periodic,  even  modulation  shown  in  Fig.  2. 

2.  When  there  is  taper  in  the  undulator  design,  the  instability  is  less  prominent  for  the  same 
reasons  as  stated  in  the  pulse  case. 


3.  Over  many  hundreds  of  passes,  the  stored  optical  wave  "sees"  many  synchrotron  oscilla¬ 
tions,  so  that  any  sideband  gain  above  threshold  gives  large  growth  from  a  small  amount  of 
noise.  The  resulting  steady-state  features  are  therefore  not  affected  by  the  details  of  the 
noise  source  employed. 

5.  The  Trapped-Particle  Instability  in  High-Current  FEL  Amplifiers 
In  the  high-current  FEL  amplifier,  y»1,  large  optical  power  can  be  produced  in  a  single 
pass  through  the  undulator.  An  induction  linac  accelerator  and  a  long  undulator  can  result  in 
values  j  =  1 04— >1 0s  [2],  The  electron  pulses  are  long,  a2  ->  so  that  the  wrapped-window  simu¬ 
lation  method  is  essential.  The  FEL  growth  rates  are  so  large  that  electrons  become  trapped  early 
in  the  undulator  and  begin  executing  synchrotron  oscillations.  There  can  be  from  several  to  many 
tens  of  synchrotron  oscillations  along  the  undulator  so  that  the  optical  field  experiences  far  fewer 
synchrotron  cycles  than  in  the  oscillator  case.  However,  even  the  limited  number  of  synchrotron 
cycles  can  result  in  significant  sideband  gain  owing  to  the  large  current  density  j. 

The  wrapped-window  simulation  in  Fig.  3  follows  the  power  and  gain  from  x-  0-»1  for  an 
FEL  with  high  current  j  =8xl04,  an  initial  field  a0=  20  starting  at  v0  =  0,  and  N  =  50.  No  electron 
phase  noise  is  present,  but  a  small  initial  sideband  field  as  =  0.01  is  introduced  at  vs  =  20jt.  Spon¬ 
taneous  emission  or  electron  shot  noise  from  the  trapped  electrons  can  lead  to  a  contribution  of 
this  size  because  j  is  large.  The  power  P  (x)  and  gain  G  (x)  are  seen  to  oscillate  with  the  synchro¬ 
tron  frequency  v5  =20jc  and  impose  a  strong  modulation  on  the  light  wave  envelope  a{z).  The 
final  power  spectrum  P  (v)  shows  that  the  sideband  has  grown  to  almost  equal  the  fundamental 
power. 
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Fig.  3.  The  trapped-particle  instability  can  create  significant  sideband  power  during  one  pass  in 
the  FEL  amplifier.  The  high-currents  used  in  amplifiers  give  much  stronger  optical  fields  at 
early  stages  of  the  undulator  that  can  cause  several  synchrotron  oscillations,  and  impose  a 
corresponding  modulation  onto  the  optical  wave  envelope.  In  this  example,  the  sideband  is 
spaced  at  vs  -  20n. 

The  trapped-particle  instability  in  FEL  amplifiers  differs  in  many  ways  from  the  oscillator 
case.  Sources  of  noise  are  much  more  important,  and  there  is  no  resonator  Q  to  consider. 

1.  Large  current  j  makes  the  trapped-particle  instability  more  likely.  Using  the  sideband 
growth  rate,  we  can  evaluate  the  current  j*  needed  to  make  the  sideband  field  as  equal  to 
the  fundamental,  a0  =  vl-  This  gives 
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The  current  j*  has  a  slow  logarithmic  dependence  on  the  initial  sideband  field  as  so  that  an 
accurate  estimate  is  not  too  crucial;  we  take  as  =  0.01.  We  also  estimate  that  there  must  be 
at  least  a  few  synchrotron  oscillations,  say  v5  =  1  6rc,  before  the  optical  wave  could  accu¬ 
rately  determine  the  synchrotron  frequency.  This  gives  a  characteristic  current  density 
j'  ~  3x1 04  where  sidebands  could  be  expected  to  significantly  alter  the  final  FEL  amplifier 
spectrum. 

2.  In  the  tapered  undulator  case,  less  current  is  trapped,  and  the  synchrotron  oscillations  have 
a  smaller  amplitude  so  that  the  trapped-particle  instability  is  observed  to  be  less  severe,  as 
in  FEL  oscillators. 

3.  Input  noise  at  the  sideband  frequency  is  important  to  the  development  of  significant  side¬ 
band  power  in  FEL  amplifiers.  Shot  noise  and  spontaneous  emission  cannot  be  eliminated, 
and  may  be  a  significant  contribution  at  large  current  densities. 

4.  For  the  large  currents  j>j*,  the  FEL  does  not  reach  steady-state  operation  even  in  strong 
optical  fields  \a  | .  The  power  continues  to  increase  and  so  does  the  synchrotron  frequency. 
This  continual  change  in  the  synchrotron  frequency  may  play  a  useful  role  in  suppressing 
the  growth  at  any  particular  sideband  frequency. 

In  conclusion,  we  find  that  the  trapped-particle  instability  is  more  difficult  to  suppress  in  the 
FEL  oscillator  than  in  the  FEL  amplifier.  The  large  number  of  synchrotron  oscillations  experienced 
by  the  light  stored  in  the  oscillator  make  the  effective  gain  over  for  many  passes  large.  But,  while 
the  trapped-particle  instability  has  already  been  observed  in  the  FEL  oscillator,  it  has  not  yet  been 
found  in  the  FEL  amplifier.  The  many  differences  between  FEL  oscillator  and  amplifier  simula¬ 
tions  implies  that  we  should  not  use  the  oscillator  case  as  a  direct  proof  that  there  will  be  an 
equally  prominent  instability  in  real  FEL  amplifiers. 
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ABSTRACT 

The  free  electron  laser  can  be  described  by  solving  the  Lorentz-Maxwell  equations 
self-consistently  in  weak  optical  fields.  The  field  evolution  is  determined  by  an  integral 
equation  that  allows  the  inclusion  of  an  arbitrary  electron  distribution  function  in  a 
simple  way.  Contour  maps  are  used  to  show  the  gain  degradation  due  to  an  electron 
beam  energy  spread  and  an  electron  beam  angular  spread.  In  the  limit  of  low  gain,  the 
gain  spectrum  is  related  to  the  spontaneous  emission  line-shape  through  successively 
higher  derivatives.  In  the  limit  of  high  gain,  it  is  shown  that  the  growth  rate  becomes 
less  susceptible  to  degradation  from  the  electron  beam  quality. 

I.  Introduction 

In  a  free-electron  laser  (FEL),  a  relativistic  electron  beam  amplifies  a  co-propagating, 
coherent  optical  wave  traveling  through  a  periodic  undulator  magnetic  field  [1].  In  the  oscillator 
configuration,  coherent  electron  bunching  develops  on  each  pass  while  resonator  mirrors  allow 
the  stored  optical  power  to  grow  over  many  passes.  In  the  amplifier  configuration,  coherent 
electron  bunches  develop  rapidly  in  the  first  part  of  the  undulator  followed  by  rapid  growth  of  the 
optical  field.  Maintaining  the  coherence  of  the  electron  bunches  over  a  significant  interaction 
length  imposes  important  restrictions  on  the  electron  beam  quality.  An  energy  or  angular  spread 
(due  to  emittance)  contributes  a  random  component  to  the  electron  motion  that  decreases  the 
coherent  bunching  in  time. 


Some  of  the  earliest  FEL  experiments  used  electron  beams  that  were  essentially 
monoenergetic  [2-5],  but  practically  all  subsequent  experiments  have  made  use  of  higher  current 
sources  with  significant  energy  spread  or  emittance.  Many  accelerators  present  a  design  trade-off 
between  high  current  and  high  beam  quality.  This  makes  it  essential  to  accurately  evaluate  the 
effects  of  beam  quality  in  present  and  future  experiments.  It  is  particularly  important  for  FELs 
designed  to  operate  at  XUV  or  X-ray  wavelengths  [6].  Several  theoretical  models  involving 
simulations  and  plasma  dispersion  relations  have  discussed  the  detrimental  effects  of  electron 
beam  quality  in  the  FEL  interaction  [7-24],  The  theory  presented  here  uses  a  convenient,  yet 
powerful,  method  of  including  an  arbitrary  electron  distribution  function  in  a  self-consistent  integral 
equation  for  the  complex  optical  field.  FEL  gain  and  the  effects  of  beam  quality  can  then  be 
calculated  analytically  or  integrated  on  a  small  computer. 

Since  the  basic  equations  solved  here  are  the  same  as  in  computer  simulations  or  the 
plasma  dispersion  methods,  specific  physical  results  have  been  shown  to  agree  with  those 
methods  when  a  direct  comparison  is  possible.  The  computer  simulations  have  proved  to  be  a 
useful  method  of  understanding  many  aspects  of  the  FEL  interaction,  but  one  of  the  most  difficult 
effects  to  accurately  characterize  is  that  of  electron  beam  quality.  Even  a  prohibitively  large 
number  of  sample  particles  is  far  short  of  the  number  in  a  real  experiment,  and  yet  introduces  a 
large  amount  of  numerical  noise  when  distributed  over  a  large  volume  of  phase-space.  To 
reproduce  some  of  the  results  shown  later  in  this  paper,  we  found  the  simulation  method  to  be 
several  hundred  to  a  thousand  times  less  efficient.  While  many  other  FEL  topics  are  most 
efficiently  studied  through  simulations,  the  detrimental  effects  of  beam  quality  are  probably  better 
handled  through  a  combination  of  analytic  and  numerical  techniques.  The  stability  analysis  used 
to  obtain  plasma  dispersion  relations  usually  calculates  the  reduced  FEL  growth  rates  due  to  poor 
beam  quality.  This  method  can  lead  to  analytical  expressions,  but  depends  upon  specific  models 
tor  the  electron  beam  distribution,  and  does  not  easily  describe  more  complicated  transient 
behavior  where  the  FEL  growth  rate  is  not  constant;  the  FEL  is  often  designed  to  operate  in  this 
regime.  In  addition,  the  exact  formulation  presented  here  works  smoothly  between  different 
regimes  of  operation  like  high  and  low  gain.  The  only  requirement  is  weak  optical  fields. 


II.  Basic  Theory 


We  solve  the  electron  Lorentz  and  optical  wave  equations  self-consistently  with  the 
assumption  of  weak  optical  fields.  The  effects  of  beam  quality  are  typically  less  important  when 
the  optical  field  strength  is  large  near  saturation,  and  the  issue  of  beam  quality  is  most  important 
in  weak  fields  where  the  accurate  evaluation  of  gain  can  determine  whether  the  FEL  is  above  or 
below  threshold. 

The  electrons  travel  through  a  periodic  undulator  with  the  field  on  the  z-axis  described  by 
Et=B  [cos0fcoz),sin(jtoz),O]  wherefi  is  the  peak  magnetic  field  amplitude.  The  undulator  field 
extends  over  a  length  L  =N’k0  with  a  number  of  periods  N,  and  wavelength  XQ  =  2TJk0.  The 
electron  velocity  in  a  perfect  helical  orbit  is  c]?  =  c  [  -  (K/y)  cos (£0  z),  ~{K/y)  sin(ifco  z),  p0  ]  where 
K  =eB  XyZranc  2,  e  is  the  electron's  charge  magnitude,  m  is  the  electron’s  mass,  c  is  the  speed  of 
light,  p0  =  (  1—  (1+/T  W)v\  and  T«c2  is  the  electron's  energy.  Imperfect  injection  due  to  poor 
beam  quality  is  more  meaningfully  introduced  after  some  further  theoretical  development.  A 
typical  undulator  uses  B  =  2kG  and  ^o  =  5cm,  so  that  K  =  1.  Since  the  electrons  are  relativistic 
(y»l),  the  transverse  excursions  ?re  small  compared  to  1$. 

The  optical  vector  potential  with  the  polarization  that  best  couples  to  the  above  trajectory  is 
X=k~'\E  |[sin\y, cosy, 0]  where  \jr  =  *z-cor+<!>,  and  X  =  2nlk  =  2rtc /co  is  the  optical  carrier 
wavelength.  The  complex  electric  field  envelope,  £(z,/)=|£(z,r)|c'w,-') ,  is  taken  to  vary  slowly 
in  z  and  r,  so  that  terms  containing  two  derivatives  in  the  wave  equation  are  small  compared  to 
terms  with  single  derivative  [25].  No  transverse  (x.y)-dependence  is  included  so  that  diffraction  is 
taken  to  be  a  small  over  the  interaction  length  L ,  and  the  electron  beam  remains  aligned  near  the 
center  of  the  optical  mode.  The  transverse  motion  above,  proportional  to  {K/y),  defines  the 
transverse  current  for  each  electron  in  the  beam.  If  the  current  density  is  uniform  over  a  sufficient 
length,  each  point  z  +  ct  in  the  optical  field  envelop  '.volves  according  to  the  slowly-varying  wave 
equation  [25] 


where  a  =  4NncKLE/yzmcz  is  the  dimensionless  optical  field  strength,  x  =  ctlL  is  the 
dimensionless  time  (Osxs  1),  j  =  8N {ncKL)zp ly3mcz  is  the  dimensionless  current  density,  p  is 
the  actual  electron  particle  density,  £  =  (k+k0)z- (at  is  the  electron  phase  in  the  combined  optical 
and  undulator  fields,  and  <...>  represents  a  normalized  average  over  all  electrons  in  the  beam 
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driving  a{x).  The  electrons  are  labeled  by  their  initial  phase-space  coordinates;  the  initial  phase  is 
C<-  =  C  (0) .  and  the  initial  phase  velocity  is  \t=d  C(0)/rf  'i  =  L[{k+k0)$0-k] .  There  are  a  large 
number  of  electrons  spread  randomly  over  each  optical  wavelength  (-107  ),  so  that  the  £,•  can  be 
accurately  taken  to  be  uniformly  spread  along  each  section  of  the  electron  beam  one  wavelength 
of  light  long,  it  can  be  easily  seen  in  (1)  that  bunching  the  electrons  near  the  relative  phase 
£  +  <{>  =  n  drives  the  optical  wave  amplitude  producing  gain,  while  bunching  near  C  +  =  n/2  drives 
the  optical  phase  <j>  without  gain.  Bunching  electrons  near  C  +  4>  =  0  results  in  negative  gain,  or 
absorption.  The  dimensionless  electron  phase  velocity  v,  has  an  initial  spread  associated  with  the 
beam  quality. 

The  electron  motion  in  the  presence  of  the  optical  wave  is  described  by  the  Lorentz  force 
equation;  dy/dt  =-{elmc)$E?.  In  the  FEL,  it  is  important  to  distinguish  between  collective 
Coulomb  forces  and  collective  high-gain  effects  [26],  Most  FELs  do  not  use  current  densities 
large  enough  for  Coulomb  forces  to  be  a  significant  effect  for  the  relativistic  electrons;  yet,  high 
gain  is  possible  and  will  be  included.  Using  the  definitions  and  assumptions  above,  the  Lorentz 
force  takes  on  the  form  of  the  pendulum  equation  [27], 


A 

dx2 


dv 

—  =  l«!cos  (C+4>) 
dx 


(2) 


The  combined  equations  (1)  and  (2)  are  valid  in  weak  or  strong  optical  fields,  for  large  or  small 
gain,  and  for  an  arbitrary  electron  distribution.  Strong  fields  near  saturation  mean  that  |  a  \  »n, 
and  weak  fields  occur  when  |a  |  «  n.  High  gain  is  achieved  when  j  »  1  and  low  gain  occurs 
when  j  <  1  [28],  Useful  FEL  configurations  display  a  wide  range  of  current  densities.  The 
electron  beam  area  is  typically  between  1mm  and  5mm,  but  the  current  ranges  from  1 A  up  to 
10kA.  Undulator  lengths  L  now  range  from  1m  to  5m,  but  will  soon  be  made  to  L  =20m  and 
beyond.  With  electron  energies  in  the  range  lOMeV  to  IGeV,  the  corresponding  values  of  j  are 
from  unity  to  more  than  5x1 04  [1],  Both  the  high-gain,  single-pass  and  the  low-gain,  oscillator 
configurations  have  important  applications. 

Equations  (1)  and  (2)  were  originally  derived  [25]  for  the  more  general  case  where  the 
electron  energy  can  change  significantly  during  a  single  pass;  in  this  case,  an  additional  factor 
H  =  (1  -\I2kN)  alters  the  wave  and  electron  equations  so  that  a  =-j <  >  and, 

C  =  v=  | a  h2cos  (C  +  b)  with  (')  =  </(  )ldx  .  The  following  work,  however,  will  be  confined  to  weak 
optical  fields  where  rj  =  1 .  An  extension  to  higher  harmonics  and  linearly  polarized  undu^tors  is 
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also  possible  without  any  change  in  form  of  (1)  and  (2),  so  that  the  general  conclusions  and 
methods  of  this  paper  are  immediately  applicable  to  a  wide  range  of  FEL  designs. 

We  now  proceed  to  solve  (1)  and  (2)  in  weak  fields,  \a  | « it,  to  obtain  an  integral  equation 
for  a[x)  incorporating  an  arbitrary  electron  distribution  function.  The  electron  phase  can  be 
expressed  as  £  =  £,•  +  v>x  +  £ (1)  where  £ (1)  is  the  first-order  perturbation  in  a .  Expanding  (1)  and 
(2)  we  have 

t 

a  (t)  =  a0  +  tijQds  <  exp[  -i  (  C  +  )  K  (1|(s )  >  .  (3) 

C  <1)('T)  =  ~J0^  jQdu[a(u)  exp(  i(  £,•  +  v,u  )  )+a*{u)  exp(-j(  £,  +  v,w  ))  J  , 

where  the  initial  optical  field  is  a(0)=  |  a  (0)  |  =a0  and  <>(0)  =  0.  We  have  made  use  of 

2k 

<  exp(  -i  £, )  >  =  JQ  d£,exp(-i£,  )/2ji  =  0  ,  since  the  initial  electrons  are  spread  uniformly  in 

phase.  The  reference  to  the  individual  electron  phases  £(1)  can  be  explicitly  removed  by 

combining  the  equations  in  (3).  Then,  we  have  an  integral  equation  governing  the  evolution  of  the 
optical  field  a  (x): 

«(t)  =  flo  +  —  \0ds\Qd<l\Qdu  <exp(-ivi(s-«))>a(«)  ,  (4) 

where  <...>  is  now  an  average  over  the  initial  electron  velocity  distribution,  and  all  reference  to  the 
electron  phases  has  been  removed.  Since  (4)  is  an  iterated  triple  integral,  it  may  be  rewritten  [29] 
as  a  double  integral, 

ii  x  * 

a(x)=a0+—jodsjodq<exp{-ivi(s-<i))>(s-q)a{q)  .  (5) 

A  normalized  electron  distribution  function  /  (vj  can  be  used  to  evaluate  the  remaining  average: 

•m  o o 

<  >  =jj*v./(v, •)(•••)  with  J_rfv,/(vi)  =  1. 

III.  Simple  Electron  Distributions 

We  begin  by  considering  two  simple  examples  with  perfect  beam  quality.  In  the  first,  we 
start  the  FEL  on  resonance  where  the  electron-optical  wave  coupling  is  largest,  /(v,)  =  8(v,). 
The  optical  wave  is  most  simply  determined  from  (4). 


(6) 


a(x)  =  a0  +  j jodsjodqjodu  a(u) 

The  integral  equation  (6)  can  also  be  written  in  a  differential  form  by  taking  successive  derivatives, 

3 

a  (x)  =  ija  (x)/2.  The  complete  solution  uses  the  form  a  (x)  =  2  a„exp(  a„x )  where  the  a„  are  the 

three  complex  roots  of  the  cubic  equation  a3  -  ij/2  =  0,  and  the  coefficients  a„  are  determined  by 
the  initial  conditions  a(0)  -a0 ,  anda(0)  =a(0)  =  0  £1].  The  solution  fora(  x )  is 

a  (X) «  —  [ exp  ( ( j/2  )'*  (i+V3 )  x/2 )  +  exp  ( ( j/2  )'*  (1-V3 )  x/2 )  +  exp  ( -  /  ( j/2  )1/3  x )  1  .  (7) 

If  the  current  density  is  small  j  ->0,orx«1  ,  we  have  the  trivial  resulta(x)  =  a0(1+yx3/12+...) 
for  an  FEL  starting  on  resonance.  There  is  no  change  in  the  optical  amplitude  | a  \  =a0+...  to 
lowest  order,  and  therefore  no  gain.  The  optical  phase  <J>(x)  increases  slowly  in  proportion  to  x3. 
The  FEL  gain  is  defined  as  G  (x)  =  (  |  a  (x)  | 2  -  a  f  )/a  f  ,  and 

G(x)  =  [2  cosh  ({j/2  )1/3  ^3  x)  +  4  cos  ( ( y/2  )1/3  3  x/2  )  cosh  ((j/2  J1*  <2x/2  )  -  6  ]  .  (8) 

In  the  high  current  limit,  j»  1  on  resonance,  the  expressions  simplify  because  one  fastest- 
growing  root  dominates  and  describes  exponential  growth  in  x.  As  seen  from  (7)  there  is  little 
change  in  the  field  during  the  bunching  time,  x  <  xfl  =  ( 2 Ij  )1/3,  that  precedes  exponential  growth. 
During  this  time,  the  electrons  move  from  their  initially  uniform  phase  distribution  to  bunch  near 
the  phase  C  +  4'  =  ^2.  As  soon  as  bunching  forms,  the  high  current  immediately  causes 
exponential  field  growth  and  high  gain.  Then, 

o(x)  =  —  exp  [  ( j/2  )1/3  V3  x  /2  ]  ,  and  G(x)  =  —  exp  [  ( j/2  )1/3  V3  x]  .  (9) 

3  9 

A  second  simple  example  is  a  high  quality  electron  beam  starting  off  resonance  at  v0.  This 
is  characterized  by  /  (v,)  =  8{  v;  -  v0 ).  The  optical  field  is  then  determined  by 

a(x)  =  a0  +  ^jodsjodqjodu  e',v,i,~u)a(u)  .  (10) 

For  low  current,  j<A,  the  optical  field  evolution  away  from  a 0  is  small  so  that  a(u)  =a0  can  can  be 
extracted  from  the  integrand  in  (10).  The  resulting  integrals  are  easily  solved  to  obtained  the 
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usual  low  current  gain  and  phase  shift  formulas  [27], 

2  -  2cos(v0t)  -  v0x  sin{v0t)  2  sin(v0x)  -  v0t  ( 1  +  cos(v0i) ) 

G(z)=j - - -  .  and  4>(x)  =  j - r -  (11) 

v0  2vo 

The  maximum  final  gain  is  G  =0.135 j  at  v0  =  2.6  and  x=1,  while  the  range  of  modes  with 
significant  gain  is  8v0  =  1  about  the  peak. 

In  order  to  obtain  the  general  solution,  use  the  substitution  b  =a  e‘v’T  in  (10).  Successive 
derivatives  then  lead  to  the  differential  form  of  (1 0),  b  -  iv0b  -  ijbl 2  =  0.  Solutions  of  the  form 

3 

b  =  £Mxp(  O-t  )  have  roots  a„  that  satisfy  the  cubic  equation  a®  -  i v0a2  -  ij/2  =  0.  In  the  limit 
«=1 

of  high  current  j  »  1 ,  the  exponential  gain  coefficient  is  reduced  by  the  factor  -  Vq  /  3^3  ( j/2  )v3 
so  that  the  gain  spectrum  is  centered  about  v0  =  0  with  a  characteristic  range  8v0  =  4.22 jv6.  In  the 
high  current  case,  the  range  of  modes  with  significant  gain  increases  slowly  as  j  increases.  We 
go  on  now  to  look  at  more  interesting  PEL  distributions  describing  less  than  perfect  electron 
beams. 

IV.  More  General  Electron  Distributions 

New  cases  of  interest  involve  more  complicated  distributions  /(v,)  describing  the  initial 
electron  beam  in  the  integral  equation  (5).  Two  electrons  starting  at  the  same  phase  £,  at  the 
beginning  of  the  undulator  (t  =  0),  but  with  different  z-velocities,  c  p0  and  c  (p0  +  Ap0),  will  drift  apart 
as  they  travel  through  the  undulator.  The  amount  of  drift  is  not  easily  predicted  without  solving  the 
full  problem,  because  electrons  can  influence  each  other  through  the  self-consistently  evolving 
optical  wave.  In  this  sense,  the  effect  of  FEL  beam  quality  is  collective.  However,  the  times  for 
the  two  electrons  to  traverse  the  undulator  are  nearly  identical  since  they  are  relativistic,  L/p0c  ~ 
Z,/(P0  +  Ap0)c  ~Lk.  An  estimate  of  their  separation  at  the  end  of  the  undulator  (t=1)  is 
Az=ApoL,  and  their  approximate  phase  difference  is  A£  =  (k  +  *0)Az  ~kAz  =kL&$0  .  If  the 
velocity  difference  Ap0  is  due  to  an  initial  energy  difference  A ymc 2  ,  we  have  Ap0  =  (l+K  2)a yiy3  , 
and  an  approximate  final  phase  separation  A£  =  4nNAy/y . 

Any  random  phase  difference  A£~  n,  or  larger,  between  electrons  in  the  beam  is  important 
to  the  FEL  operation,  because  bunching  on  the  optical  wavelength  scale  is  diminished 
significantly.  At  the  end  of  the  undulator,  the  final  phase  difference  is  roughly  estimated  by 
A^  =  Av(  for  each  electron.  From  the  definition  of  the  electron  phase  velocity  v;  ,  we  see  that  a 
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small  change  in  the  initial  electron  energy  Ay  me 2  corresponds  to  a  change  in  the  initial  phase 
velocity,  AVj  =  4jiN Ay  /  y  .  y  »  1 .  A  distribution  of  initial  electron  energies  from  an  accelerator  or 
storage  ring  is  often  accurately  represented  by  the  normal  distribution  function  so  that  we  can  take 


/(v,)  = 


exp(-(v,-v0)2/2CT2) 

V2ita 


(12) 


where  a  is  the  standard  deviation  of  V;  away  from  the  peak  phase  velocity  v0.  If  Ay  me 2  is  taken  to 
be  the  standard  deviation  of  the  electron  energy  away  from  y me2,  then  c~4nNAyly  .  Two 
electrons  starting  at  the  same  phase  ,  but  with  an  energy  difference  Ay=  yl  4N  will  drift  apart  by 
roughly  half  of  one  optical  wavelength  at  the  end  of  the  undulator.  A  random  spread  of  width  a = n 
causes  a  random  phase  spread  of  approximately  AC  =  n  at  the  end  of  the  undulator  and  impairs 
bunching.  Inserting  (12)  into  (5)  gives 


a{T)  =  a0+^f  dsf  dq  e^^'12  (s-q)a(q)  . 


The  gaussian  factor  in  the  integrand  decreases  the  coupling  current  j  as  t  increases,  and 
describes  the  degradation  of  bunching  due  to  the  spread  in  electron  phase  velocities.  The 
complicated  self-consistent  evolution  of  the  electron  beam  distribution  and  the  optical  field  are 
described  exactly  in  (13),  but  before  evaluation,  we  can  generalize  its  form  further. 

An  angular  spread  is  also  possible  due  to  the  finite  emittance  of  an  electron  beam.  An 
electron  of  energy  y me2  entering  the  undulator  with  a  small  injection  angle  0,  has  a  reduced  z- 
velocity,  p0  -» p0  cos  0,  =  p0  ( 1  —  0 2  /  2 )  .  The  resulting  z -velocity  change  is  Ap0  =  -  G,2  /  2  , 
reducing  the  initial  phase  velocity  by  Av,  =  -  2nNy  202/(1+/:  2 ) .  A  gaussian  distribution  of  angles 
about  the  z -axis  with  standard  deviation  AS,  gives  the  exponential  distribution  function 


exp  ( -  ( v0  -  v, )  /  cre ) 

/  (v,)  = -  for  v,  <v0  ,  and  /(v,)  =  0  for  v,  >v0  ,  (14) 


where  a9=  47iA/y2A0,  2/(1+/f  2),  and  v0  is  the  phase  velocity  for  electrons  entering  on-axis.  The 
distribution  function  (14)  is  sharply  peaked  at  v<  =v0  where  electrons  enter  on-axis,  and  decays 
exponentially  for  v,  <  v0  because  the  injection  at  any  angle  8,  can  only  decrease  the  electron’s 
longitudinal  velocity  and  its  phase  velocity  V;.  If  each  element  of  the  energy  distribution  (12)  is 
given  an  angular  spread  according  to  (14),  then  the  resulting  integral  equation  for  the  optical  field 
becomes 
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(15) 


The  transverse  motion  of  electrons  injected  at  an  angle  is  sometimes  confined  by  either  the 
natural  off-axis  undulator  fields  or  external  focusing  elements.  The  focusing  forces  result  in 
transverse  betatron  oscillations  about  the  undulator  axis.  When  the  electron  beam  is  injected  to 
match  the  natural  focusing  properties  of  an  undulator,  the  number  of  betatron  oscillations  along 
the  undulator  is  N$~NkH2 y.  In  the  limit  of  large  y  and/or  small  K,  the  angular  spread  of 
electrons  can  be  important,  oe  >  1,  while  the  transverse  focusing  can  be  made  negligible,  <  i  . 
In  this  limit,  the  integral  equation  (15)  applies. 

The  complex  optical  field  a  (x )  now  depends  on  an  input  energy  spread  characterized  by  cr , 
and  an  input  angular  spread  characterized  by  cr8  .  Other  types  of  distribution  functions  can  be 
added  in  a  similar  way.  If  (5)  is  solved  numerically,  even  experimental  distribution  functions 
peculiar  to  a  given  accelerator  or  transport  system  can  be  added.  The  general  result  (5)  and  the 
specific  example  (15)  are  important  results  of  this  paper.  They  provide  analytic  expressions 
describing  FEL  performance  with  an  arbitrary  electron  distribution  function. 


V.  Low-Current  FELs 

One  of  the  cases  of  general  interest  is  the  low  current  FEL  oscillator.  Radiation  energy  is 
stored  in  an  optical  resonator,  and  repeatedly  driven  by  successive  electron  pulses  from  an 
accelerator  like  a  linac  or  storage  ring.  An  important  issue  for  the  oscillator  is  the  detrimental 
effect  of  the  electron  energy  and  angular  distributions  when  the  oscillator  is  starting  from  weak 
optical  fields.  In  the  low-current  case,  we  can  simplify  (15)  by  taking  a(q)  ~  a0  in  the  integrand  of 
the  integral  equation,  and  neglecting  higher-order  terms  in  j  .  Without  a(q)  in  the  integrand,  the 
integral  can  be  further  simplified  by  noticing  that  iqe~‘y’<’  =  -  3Vi  Then, 
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(16) 


Direct  integration  of  (16)  is  possible,  but  the  result  is  a  complicated  expression  containing  many 
error-functions  (30). 

We  can  alter  the  form  of  (16),  however,  to  obtain  some  important  physical  interpretations. 
The  factors  e~a"12  and  ( 1  - iaBq  )~'  can  be  interpreted  as  power  series  expansions  in  q 
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multiplying  the  factor  e  *  ;  then  {-iq )  e  -» 3;§  e  .  Now,  write 


a(x)-a  o 


exp(cr232/2) 
1  +  oe  3V> 


The  double  integral  is  simply  vq2  ( 1  -  /v0t  -  e~<v,x ),  so  that  the  complex  integrations  in  (16)  have 
been  replaced  with  a  power  series  expansion  to  all  orders  in  a  and  oe-  To  first  order  in  a  and  ct0, 
an  explicit  expression  for  a(x)  is  easily  obtained  from  (17).  This  is  a  useful  limiting  case  since  a 
low-current  FEL  system  would  not  typically  use  a  low  quality  electron  beam  ( large  o  or  ae )  and 
remain  above  threshold.  From  (17),  the  low-current  FEL  gain  at  the  end  of  the  undulator  is 

exp(o232  /2)  f  sin2(  Vn/2 ) 


We  recognize  the  factor  in  brackets  [...]  as  the  FEL  spontaneous  emission  line-shape  for  an 
electron  in  a  perfect  trajectory  through  the  undulator.  It  has  been  known  for  some  time  that  the 
gain  is  fundamentally  related  to  the  derivative  of  the  spontaneous  emission  line-shape  [31].  The 
new  feature  presented  in  (18)  is  to  express  how  the  electron  beam  energy  and  angular  spreads 
affect  that  relationship  through  successively  higher  derivatives. 

With  the  physical  interpretation  of  the  line-shape  factor  [...]  ,  we  can  substitute  alternate 
forms.  One  convenient  choice  is  [...]-» exp  (-v0z/47t)  which  approximately  reproduces  the 
correct  features  of  the  simple  FEL  gain  spectrum,  G  =  ( j v0 14n )  exp  ( -vf/4n ).  The  successive 
derivatives  evaluating  the  effects  of  beam  quality  lead  to  more  compact  expressions,  and  illustrate 
how  (18)  can  be  used  in  practical  situations.  Even  an  experimental  line-shape  could  be  used  in 
(18). 

While  the  analytic  results  presented  have  their  merit,  the  complete  integral  (15)  is  easy  to 
integrate  on  a  small  computer.  The  values  needed  for  the  contour  plots  of  this  paper  were 
evaluated  in  this  way.  Figure  1  shows  a  combined  intensity  and  contour  plot  of  ln(1+G(a,v0)) 
where  the  final  gain  at  the  end  of  the  FEL  undulator  is  G  =  (a*(1)a(1)-ag  )/a0  .  In  Fig.  1 
a9  =  0  ,  so  that  gain  degradation  is  only  due  to  an  energy  spread  with  no  angular  spread.  The 
current  density  is  7  =  1.  and  gives  low  gain  so  that  ln(1+G)  =  G  .  The  brightest  points  (white)  on 
the  (a,v0)-surface  indicate  peak  gain  G  =  0.13/,  while  the  darkest  points  (black)  indicate  maximum 
absorption  G  =  -0.13./.  Zero  gain  is  indicated  by  the  intermediate  grey  shown  in  the  scale  at  the 
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top.  Specific  contours  of  constant  gain,  ln(1+G)  =  ±0.06,±0.08,±0.1,  and  ±0.12,  are 
superimposed  on  the  intensity  plot.  The  gain  surface  is  approximately  antisymmetric  about  v0=  0, 
and  in  the  limit  j  ->  0,  the  gain  G  (a,v0)  becomes  exactly  antisymmetric.  The  characteristic 
amount  of  spread  required  to  decrease  the  gain  is  seen  from  Fig.  1  and  (18)  to  be  a  1.  Note 
also  that  as  the  spread  a  increases,  the  phase  velocity  for  peak  gain  v0*  =  2.6  increases  slightly. 
Peak  absorption  occurs  at  -  v0*  and  slightly  decreases  with  increasing  o. 


Fig.  1.  Intensity  and  contour  plot  of  ln(1+G(o,v0))  with  j  =  1  and  ae  =  0.  The  weak-field  gain 
degradation  in  this  low-current  FEL  is  due  to  an  electron  beam  energy  spread  with  a  normal 
distribution  function. 

Figure  2  shows  a  combined  intensity  and  contour  plot  of  ln(1+C(cr0,vo))  evaluated  by  (15) 
with  j  =  1  and  a = 0.  The  gain  degradation  here  is  due  to  a  monoenergetic  electron  beam  entering 
the  undulator  with  an  angular  spread  described  by  c9.  The  grey  scale  and  contours  of  gain  are 
the  same  as  in  Fig.  1.  Unlike  Fig.  1,  the  absorption  contours  (white)  have  a  much  different  shape 
than  the  gain  contours  (black).  Since  the  distribution  function  / (v,)  due  to  an  angular  spread  is 
skewed,  there  is  no  reason  to  expect  the  antisymmetric  properties  of  G(oe-O,v0)  to  be 
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maintained  at  cr0  >  0  .  As  seen  in  Fig.  1  and  (18),  the  characteristic  value  tor  the  degradation  of 
gain  is  c0*=1.  The  phase  velocity  for  peak  gain,  and  peak  absorption,  both  increase  with 
increasing  o0  roughly  as  v0*  =c0.  Note  that  the  general  features  of  Figs.  1  and  2  are  quite 
different  owing  to  the  different  forms  of  the  electron  distributions.  This  emphasizes  the  importance 
of  the  shape  of  the  electron  distribution  in  evaluating  gain  degradation  in  FELs,  and  the  need  for 
an  accurate,  flexible  theory  as  presented  here. 


XXX*X*  FEL  Gain  Surface  ****** 


Fig.  2.  Intensity  and  contour  plot  of  ln(1+G(cT9,v0))  with  j  =  1  and  <r  =  0.  The  weak-field  gain 
degradation  in  this  low-current  FEL  is  due  to  an  electron  beam  angular  spread  with  a  normal 
distribution;  the  resultant  phase-velocity  distribution  is  the  exponential  distribution  function. 

VI.  High-Current  FELs 

The  integral  representation  of  the  optical  field  in  (15)  is  also  valid  for  high-current  FELs 
where  j  »  1.  In  this  case,  a(z)  acquires  a  non-linear  dependence  on  j  ( recall  expression  (7)  ) 
and  cannot  be  removed  from  the  integrand  of  (15).  To  proceed  analytically,  it  is  convenient  to 
remove  one  integral  from  (15)  by  taking  the  t-derivative  of  both  sides,  and  use  the  form  a  =a0cax 


(or  the  optical  field.  Since  j  »  1,  assume  that  a  has  some  large  real  part,  even  though  the 
exponential  growth  may  be  somewhat  diminished  by  the  presence  of  a  and  ag.  The  magnitude  a0 
cancels  on  both  sides,  and  a  change  of  variables  gives  the  form 


«  U  t  .  -m  e  -iv.sla  , 

a  ~ — J  ds  s  e  - ~e  .  (19) 

2  1-iaes/a 

The  upper  integration  limit  in  (19)  has  been  extended  to  infinity  because  the  integrand  containing 
the  factor  e~*  becomes  negligible  for  large  s . 

Eqn.  (19)  describes  several  properties  of  high-gain  FELs  without  integration.  If  a,  v0,  and  a6 
all  -*  0  ,  then  a  has  the  same  roots  found  in  (7).  If  the  current  density  j  -> »  so  that  a  real  part  of 
a  -» <*>,  then  we  obtain  the  same  limit,  since  a,  v0,  and  ae  all  appear  divided  by  a  in  (19).  Unlike 
the  low-current  FEL,  the  importance  of  beam  quality  in  a  high-current  FEL  depends  on  the  current 
density  j  .  This  feature  has  been  seen  in  FEL  experiments  and  simulations,  but  is  now  expressed 
analytically.  The  importance  of  beam  quality  can  be  made  more  quantitative  by  iterating  (19). 
Estimating  the  real  pari  of  the  fast-growing  root  as  a*  =  ( j/2  )1/3V3/2,  the  integrand  of  (19)  is  only 
significantly  modified  when  a*  =  (j!2)m  ^3/2  or  when  <sg*  ~{jl2)v3^3/2.  In  the  high-current 
FEL,  the  characteristic  values  of  beam  quality,  a  *  and  cr8* ,  are  not  equal,  and  increase  with  the 
current  density  j.  These  expressions  should  be  helpful  in  designing  high-gain  experiments  where 
there  is  a  trade-off  between  beam  quality  and  beam  current. 

Figure  3  shows  a  combined  intensity  and  contour  plot  of  In  (1+G  (cr, v0))  for  moderately  high 
current ;  =  1 00  and  ae  =  0.  The  points  at  the  peak  gain  In  (1+G )  =  4.3  are  indicated  by  white  on  the 
(ct,v0) -surface;  black  indicates  zero  gain.  Contours  of  constant  gain,  In  (1+G )  =  2.0,  2.5,  3.0,  3.5, 
and  4.0  are  superimposed  on  the  intensity  plot.  For  the  high  quality  electron  beam,  small  a,  gain 
is  confined  to  a  region  near  resonance,  but  extends  to  a  broader  range  in  v0  than  in  the  low- 
current  cases  of  Figs.  1  or  2.  This  agrees  with  the  discussion  below  (11),  and  gives  the  range  of 
optical  wavelengths  over  which  there  is  significant  gain  Sv0~4jv6~7  .  To  find  the  range  of 
wavelengths,  use  &X  fk ~ Av0 /2nN  about  the  resonant  wavelength  X  =  X0(1+A'  z)/2y2  .  The 
maximum  available  gain  decreases  significantly  as  a  -+  a  *  ~  4.5  as  predicted  in  the  previous 
paragraph,  and  the  phase  velocity  for  peak  gain  roughly  follows  v0*  ~a  . 
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Fig.  3.  Intensity  and  contour  plot  of  In  (1+G  (o,v0))  with  j  =  1 00  and  ctb  =  0.  The  weak-field  gain 
degradation  in  this  FEL  is  due  to  an  electron  beam  with  a  normal  distribution  in  energy. 
Figure  4  shows  the  plot  of  ln(l+G(cr0,vo))  with  a=  0  so  that  the  gain  degradation  is  caused 
by  an  angular  spread  in  the  electron  beam.  The  contours  of  constant  gain  differ  from  Fig.  3 
because  of  the  new  shape  of  the  electron  distribution  function.  As  cre  increases,  there  is  a  slower 
decrease  in  gain  because  <re*  >  o  *  as  found  above.  When  expressed  in  dimensionless  form,  an 
angular  spread  is  better  tolerated  in  an  FEL  than  is  an  energy  spread.  The  points  of  peak  gain 
increase  with  increasing  cre  similar  to  Fig.  3. 
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Fig.  4.  Intensity  and  contour  plot  of  In  (1+G  (ae,v0))  with  y  =  100  and  cr  =  0.  The  weak-field  gain 
degradation  in  this  FEL  is  due  to  an  electron  beam  with  an  angular  spread  producing  an 
exponential  distribution  in  phase  velocities. 

Figure  5  shows  the  combined  intensity  and  contour  plot  of  In  (1+G  (a, v0))  for  high  current 
j  =  104  with  ae  =  0.  There  are  no  negative  gain  regions,  and  the  available  peak  gain  is  much 
larger  than  for  the  lower  current.  For  a  =  0,  the  position  of  peak  gain  is  essentially  at  resonance 
vo=0,  but  again  increases  roughly  as  v0*  ~a  while  beam  quality  diminishes.  The  width  of  the 
gain  spectrum  at  a  =  0  is  wider  than  the  lower  current  case,  and  agrees  well  with  5v0  =  4 jv6  ~  12  . 

The  contours  of  constant  gain,  In  (1+G)  =  14 . 24,  show  that  the  range  of  wavelengths  for  gain 

becomes  narrower  as  o  increases,  and  the  maximum  available  gain  decreases  significantly  as 
a— >cr*  =20. 
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Fig.  5.  intensity  and  contour  plot  of  In (1+G (a,v0))  with  j  =  104  and  cre=  0.  The  weak-field  gain 
degradation  in  this  high-gain  FEL  is  due  to  an  electron  beam  with  a  normal  distribution  in 
energy. 

Figure  6  shows  the  plot  of  ln(1+G(ae,v0))  for  high  current  j  =  104  with  a  =  0.  Again,  the 
contours  of  constant  gain.  In  (1+G)  =  14, ....  24,  are  distinct  from  Fig.  5  showing  the  importance  of 
the  electron  beam  distribution  function  even  at  high  gain.  As  ce  ->  cxe* ,  the  gain  decreases 
significantly,  but  again  the  angular  spread  is  seen  to  be  less  harmful  than  an  energy  spread. 
Unlike  Fig.  5.  the  position  of  peak  gain  stays  closer  to  resonance  as  a0  increases. 
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XXXXX*  FEL  Gain  Surface  XXXXXX 


Fig.  6.  Intensity  and  contour  plot  of  ln(1+C(a9,v0))  with  j  =  104  and  c?  =  0.  The  weak-field  gain 
degradation  in  this  high-gain  FEL  is  due  to  an  electron  beam  with  an  angular  spread 
producing  an  exponential  distribution  in  phase  velocities. 
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